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£^ ■ An X-ray flux limited sample of nearby clusters of galaxies has been defined 

C • . for observations of the Sunyaev-Zeldovich effect (SZE) to be carried out on the 

Owens Valley 5. 5- meter telescope at 32 GHz. The X-ray sample selection minimizes 
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the systematic errors introduced by cluster elongation in the determination of Hq. 
Due to their proximity, these clusters are well-studied in the X-ray wave-bands. 



o 

The measurement of the SZE in three of these clusters is reported in this paper: 



AT = -375 ± 28 //K (A478), -437 ± 25 //K (A2142) and -243 ± 29 //K (A2256). 
These values have been corrected for radio source contamination, but have not been 
corrected for the beam dilution and switching (which are model-dependent). There 
is an additional overall calibration uncertainty of 7%. If the temperature profile of 
the clusters is known, then the SZE provides a direct probe of the baryonic mass 
in the hot ionized phase of the medium. We find surface baryonic mass densities of 
(7.5 ± 2.5) x 10 13 M©Mpc~ 2 within the T.35 FWHM Gaussian beam of the 5.5-m 
telescope projected on the cluster centers. For A2142, A2256, and the Coma cluster 
previously observed by Herbig et al. (1995), we find a consistent value for the 
ratio of the SZE determined baryonic mass to the gravitational binding mass of 
M sze /M tot = 0.061 ± 0.011 h" 1 . Note that this is a lower limit on the total baryon 
fraction, as there may be significant contributions from other baryons. Comparison 
with the standard big-bang nucleosynthesis prediction f^/i 2 = 0.013 ± 0.02 gives a 
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value for the cosmological density parameter of QqIi <J 0.21 ± 0.05, assuming our limit 
on the baryon fraction in clusters applies to the Universe as a whole. This density is 
in agreement with independently determined values from large-scale structure studies. 
Recent values for f2#/i 2 based upon deuterium abundances are outside the previously 
accepted range, and in combination with our data lead to significantly higher or lower 
J)o- Finally, we present preliminary determinations of the Hubble constant using X-ray 
models gleaned from the literature. The data from the three clusters, along with the 
results previously obtained using the same telescope on the Coma cluster, yield a 
sample average value Hq = 54 ± 14kms _1 Mpc -1 . We discuss the uncertainties in 
these results and future prospects for this method. 

Subject headings: cosmic microwave background — distance scale — cosmology: 
observations — galaxies: clusters — intergalactic medium — dark matter 

1. Introduction 

The scattering of the cosmic microwave background radiation by hot gas in clusters of 
galaxies, known as the Sunyaev-Zeldovich Effect (SZE), has been recognized for over two decades 
as a potentially important tool for cosmological and astrophysical studies (Sunyaev & Zeldovich 
1980). In the SZE, inverse Compton scattering boosts the microwave background photons to 
higher frequencies, up-shifting and distorting the Planck blackbody spectrum. At low frequencies 
where the spectrum rises with frequency, this reduces the intensity of the CMB at a given 
frequency. At higher frequencies where the spectrum falls with frequency, the SZE increases the 
intensity. Relative to the background, clusters look dark at low frequencies and bright at high 
frequencies. 

A particularly important application of SZE observations is the determination of the Hubble 
constant, Hq. Classical determinations of Hq rely upon the cosmic distance ladder (e.g., Tully 
1988) and are therefore subject to the uncertainties inherent in each step of the ladder. For this 
reason "direct" determinations of Hq would be very important if the systematic errors could 
be understood and allowed for. Examples of such direct methods are expanding supernova 
photospheres, gravitational lenses and the SZE. For recent determinations of Hq using classical 
methods, see Kennicutt, Freedman & Mould (1995) and references therein, as well as Sandage 
(1996), and Mould et al. (1995). 

The high temperature ionized cluster medium produces both Compton scattering (SZE) 
and thermal bremsstrahlung (X-ray) emission that depend upon different powers of the electron 
density and temperature, n e T e for the SZE, and approximately n 2 e T^ 2 for the bremsstrahlung 
component of the X-ray. The integral equations for the observed SZE and X-ray brightness can, 
given a suitable model for the density profile and knowledge of the electron gas temperature, be 
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solved for central density and the linear core radius — when combined with the observed angular 
core radius this yields a value for the angular diameter distance D a . For clusters at low redshift, 
this yields the Hubble constant Hq. For a set of clusters covering a wide range of redshifts, the 
determination of D a (z) as a function of redshift z can constrain or determine other cosmological 
parameters such as the density parameter fio and cosmological constant Ao. 

The gas in a massive galaxy cluster has a temperature of roughly 10 s K and a central density 
in excess of 10~ 3 cm~ 3 . This leads to an expected microwave decrement along the line of sight 
through the cluster center in the range O.lmK+lmK. Despite the large size of the SZE relative 
to intrinsic anisotropy signals 50 /iK), observations of the effect have proven difficult and have 
been plagued by systematic errors. However, a number of experiments have now produced reliable 
SZE measurements. Using the 40-m telescope at Owens Valley Radio Observatory (OVRO) at 
20 GHz, Birkinshaw (1990) found AT RJ of -444 ± 65 fiK, -301 ± 49 fiK, and -354 ± 43 fiK for 
the clusters 0016+16, Abell 665, and Abell 2218, respectively. Herbig et al. (1995) measured 
a decrement of —308 ± 51 /jK in the Coma cluster using the OVRO 5.5m telescope. Recent 
interferometric measurements of the SZE in A2218 (Jones et al. 1993) and A773 (Grainge et al. 
1994) using the Ryle telescope at 15 GHz have been reported. Interferometric measurements of 
0016+16 and A773 have also been made by Carlstrom, Joy & Grego (1996) with the OVRO 
Millimeter Array outfitted with 32 GHz receivers. The SZE in the cluster A2163 has been 
measured using a bolometer array on the CSO (Wilbanks et al. 1994) — this measurement of 
the effect as an increment on the high-frequency side of the blackbody peak is the first to place 
interesting limits on the peculiar motion of massive clusters. 

These observations have all been made on fairly distant clusters. This has been dictated by 
the use of large single dishes or large interferometers. These instruments have small switching 
angles or small fields-of-view that resolve out the structure on the angular scales corresponding to 
nearby clusters. In the determination of Hq, the models of the X-ray emitting gas in the cluster 
potential are just as important as the SZE measurements. Accurate models are available only for 
the nearby clusters. Also, the control of biases to the Hq determination, such as cluster elongation 
or substructure, requires the systematic study of a complete sample of carefully selected clusters. 
At the present, the best samples for this study are drawn from nearby (z < 0.1) X-ray selected 
cluster catalogues. 

The newly-completed 5.5-m telescope at OVRO is an ideal instrument for SZE measurements 
in a sample of nearby clusters. The telescope operates at 32 GHz and has a sensitive wide- 
bandwidth HEMT receiver. The primary beam is 7'.35 (FWHM) and the dual-horn switching 
angle is 22f 16. This telescope has the right beamwidth and sensitivity to measure the SZE in 
clusters with angular core radii in the range l' to 22' with reasonable efficiency. The Herbig et al. 
(1995) measurements of the Coma cluster (z = 0.023) recovered approximately 61% of the central 
decrement. 

We report here the first SZE observations undertaken with the specific aim of estimating and, 
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where possible, eliminating the major known sources of systematic error in the determination of 
Hq by means of SZE observations of a sample of nearby clusters. 

2. The X-ray Cluster Sample 

There are three potential problems in SZE determinations of Hq: elongation, density 
substructure, and non-isothermality of the cluster gas. If a cluster is elongated along the line of 
sight, then Hq will be underestimated, while if the cluster is elongated in the plane of the sky, 
then Hq will be overestimated. Clumping of the intracluster medium also causes a bias, due to the 
different dependences of the X-ray emission and SZE upon the electron density and temperature. 
In this case, Hq is overestimated by the factor (n 2 )/(n e ) 2 , a quantity which is always greater than 
unity. The measurement of Hq is proportional to T^ 2 , so errors in temperature or temperature 
gradients will be a problem. Detailed discussion of these problems have been presented in the 
literature (e.g. Birkinshaw et al. 1991). 

The uncertainties due to substructure and isothermality can only be resolved through high 
resolution X-ray observations and detailed modeling. For this to be possible, nearby clusters must 
be chosen — new measurements by ASCA and ROSAT of these objects are revolutionizing our 
understanding of cluster astrophysics. However, if clusters are prolate or oblate ellipsoids, axial 
ratios of 0.5 imply an uncertainty of up to a factor of 2 in the line-of-sight distance through the 
cluster. This uncertainty will average out in determinations of Hq if the SZE is measured in a 
large enough orientation unbiased sample of clusters. A sample chosen by central X-ray surface 
brightness will be systematically biased towards clusters that have long axes along the line of 
sight. Such a bias arises because these clusters will have a greater central brightness than clusters 
at the same distance with their short axes along the line of sight. To guard against this selection 
effect it is necessary to choose a complete flux-limited parent sample while staying sufficiently 
above the sample flux limit to minimize the biasing effect upon the selection. 

We selected the X-ray flux- limited sample of Edge et al. (1990) for our parent sample. This 
sample is complete for fluxes in the 2-10 keV band / x > 3.1 x 10~ n erg cm -2 s -1 , and is 70%-90% 
complete at f x > 1.7 x 10~ n erg cm -2 s _1 . This sample was derived by Edge et al. from HEAO-1 
and Ariel V surveys and cross-checked with Einstein and EXOSAT observations. 

We restrict our observations to the higher flux limit f x > 3.1 x 10~ n erg cm~ 2 s , 
galactic latitude \b\ > 20°, and 5 > —23°. A high-luminosity sub-sample was selected with 
L x > 1.25 x 10 44 /i" 2 ergs _1 (2-10 keV), where Hq = 100 /ikrns" 1 Mpc" 1 . The 11 clusters in this 
subsample are listed in Table |]. The core radii for these clusters fall within the optimum size 
range for the 5.5-m telescope. Also, given the X-ray parameters for these clusters, we expect SZE 
decrements in the range 250 /iK - 750 fiK, which are easily observable with our system sensitivity. 
The redshifts and thus the luminosities listed in Table |l] were taken from the Edge et al. (1990) 
paper. 
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The Edge et al. sample for f x > 3.1 x 10 -11 erg cm -2 s -1 is complete for z < 0.066 
at L x > 1.5 x 10 44 h~ 2 ergs -1 . The lowest luminosity cluster in our sample is Coma 
[L x = 1.85 x 10 44 /i _2 ergs _1 ) — at this luminosity, the sample would be complete for z < 0.073. 
The most distant cluster in our sample is A2142 {z = 0.09); at this redshift the sample would 
be complete for L x > 2.82 x 10 44 h~ 2 ergs -1 . In order to be able to observe a reasonably large 
number of clusters, we choose to use the entire sample, though it is flux-limited only, and 
not volume-limited out to z = 0.1. This means, therefore, that we are prone at the lowest 
luminosity levels to selection effects such as elongation, and some care will have to be taken in the 
interpretation of the statistical results. 

From the parent sample, we have selected for our first observations clusters that are free of 
strong radio sources which would contaminate our SZE measurements at 32 GHz. The clusters 
A478 and A2142, massive clusters with well-measured X-ray profiles, were ideal targets, with large 
expected SZE decrements. The cluster A2256 is optimal because of its high declination. The 
Coma cluster was previously observed with this instrument by Herbig et al. (1995), and has been 
included in the sample. 



3. Observations 

Although SZE signals are roughly an order of magnitude larger than the limits that have 
been placed on intrinsic anisotropy signals, the removal of systematic effects from SZE data still 
requires great care. The SZE measurements must be very sensitive, with uncertainties less than 
10% in order to contribute less than 20% to the uncertainty in the estimate of Hq. Most clusters 
are tracked over large angular distances across the sky. The variation of azimuth and zenith angle 
can introduce significant ground spillover effects, which can be minimized in intrinsic anisotropy 
observations by observing only fields near the North Celestial Pole. Our observations of clusters 
use a three-level differencing technique to remove systematic effects from ground spillover, the 
atmosphere, and the receiver. 

The observations reported here were made between 15 July 1993 and 6 March 1994. The 5.5-m 
telescope at OVRO was used, employing a receiver operating at a center frequency of 32 GHz, with 
a primary beam FWHM of 7'.35. The receiver noise level for a measurement was approximately 
0.9mK in one second of total integration time. The noise level including atmosphere and ground 



pick-up is about 50% higher than this. The receiver is discussed in § 3.1 and the sensitivity and 
noise performance in 



This section contains detailed descriptions of the observing, calibration, data editing, filtering 
and analysis procedures. Much of the discussion in may be skipped by the casual reader, though it 



is suggested that the reader look at the discussion of the referencing in § 3.2 and the data analysis 



in § 3.5 before proceeding to the presentation of the results for the cluster observations in 
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3.1. First and Second Differencing 

The receiver has two horns: ANT (antenna) and REF (reference). These have nearly identical 
Gaussian on-sky response patterns with 7'35 FWHM and a separation of 22fl6. The rms pointing 
accuracy of the 5.5-m telescope was found to be better than 0'.5 from observations of bright 
calibration sources. The receiver switches between the two beams every 1 ms (/ = 500 Hz) by 
means of a pulse latched switch. The output of the switch then passes through an amplifier chain 
before detection. Our observations employed as the first stage a TRW InSb HEMT amplifier with 
a bandwidth of 5.7 GHz centered at 32 GHz. After detection, the power measurements taken in 
each switch orientation are accumulated and the mean difference and the standard deviation of 
the mean difference are computed for each I s fundamental integration period. This first level of 
switching removes the offsets and signal common to both beams, and cancels the offsets due to 
low-frequency gain fluctuations in the amplifier chain. 

Another level of switching is required to remove gradients. The telescope is moved in azimuth 
by the switching angle 22fl6 and a second difference is formed. The period for this switching is 
20 s to 50 s . The second differencing removes gradients in the atmospheric emission that are stable 
on these timescales. 

The two 5.5-m beams of the telescope cease to intersect at a height of 853 meters for the angle 
of 22' 16. For cloudlets moving through the beams, a characteristic speed of 1 m/s transports the 
fluctuations through the beams on a timescales of 5 s at this height. The 500 Hz switching "freezes" 
these fluctuations, which then cancel if they move through both the ANT and REF beams. We 
expect these fluctuations to cause an increased noise level in the second difference measurements. 

The basic measurement is the "FLUX" (see Readhead et al. 1989 or Myers et al. 1993 for a 
description of this procedure). In a FLUX measurement, the telescope first moves (slew time r s ) 
to a position where the REF beam is pointed at the cluster (position "ON") and the ANT beam is 
22'16 away in azimuth (position "Rl"), where it spends time r« integrating. The single-differenced 
power recorded in this position is denoted "A". The telescope is then moved (r s again) so that 
the ANT beam points at the cluster ("ON" again) and the REF beam is displaced by 22'16 in 
the opposite direction (position "R2" ) from the initial ANT position. Two consecutive integration 
cycles are performed in the new position, yielding power measurements "B" and "C" . Finally, the 
telescope is moved back to the original position for an additional integration, denoted "D" . If a 
series of FLUX measurements is to be performed consecutively on the same target, the subsequent 
"A" integrations are started from the same position as the previous "D" integration, without a 
telescope slew. Therefore, each FLUX beyond the first in a series takes a total time 

T f = A Ti + 2T s , (1) 

with the first FLUX in the series taking r s longer for the first slew. Typically, t s = 12 s and 
Tj = 20 s , so a FLUX lasts tj = 104 s . The actual slew time t s depends upon the zenith angle of 
the observation, increasing from 12 s for zenith angles greater than 30°, to more than 30 s close to 



-7- 



the zenith. 

The individual integrations are normalized and calibrated such that they measure the 
antenna temperature differences {Tant — Tref)- We therefore construct an individual FLUX 
measurement 

FLUX= -(-A + B + C- D). (2) 

The telescope also records a standard deviation (SD) measurement for each FLUX by measuring 
the variance of the I s difference measurements during the integrations. This procedure subtracts 
the average of the signals in two reference positions on either side of the cluster from the signal on 
the cluster. In terms of antenna temperature on the sky, 

FLUX = T ON -\ (Tm + T R2 ) . (3) 

We can also rearrange the individual integrations A-D, which are recorded separately as part 
of the data log, to form other quantities of interest. For example, we define the switched difference 

SW1 = ^(-A-B + C + D). (4) 

Because the differences (D — A) and (C — B) appear in the SWl, and are the integrations at the 
same position, the signal from the sky in the far-field is cancelled, leaving only time variations of 
the atmosphere (possibly drifting through the slightly divergent beams), and spurious signals not 
common to the integration pairs. There is another combination of the integrations that cancels 
out the far-field signal, 

SW2 = -(-A + B -C + D), (5) 

which also acts like a time-difference filter, but with twice the frequency of the SWl. The more 
moderate filtering provided by the SWl has proved to be a useful diagnostic for periods of bad 
data, without any bias due to strong celestial sources in the beams. Note also that, because of the 
form of (|j) and (||), the SD from the FLUX applies to these quantities also. However, unlike the 
FLUX, the SWl and SW2 are only single-differences in time, and tend to have correspondingly 
larger scatter about the mean values than the FLUXes. These quantities are useful for throwing 
out grossly discrepant measurements but do not provide a stringent filter for more subtle problems. 

As the telescope tracks the cluster across the sky, the reference beams trace circular arcs 
around the cluster. The reference beams are always separated from the ON beam purely in 
azimuth. The position of the reference beams for a given FLUX is recorded in the form of the 
parallactic angle (ip), defined as the angle between the direction of the North Celestial Pole and 
the direction of the zenith. The parallactic angle is 

, cos A sin if 

tan^>=^- -z r-r-j 77, (6) 

sin A cos o — cos A sin o cos ti 

where A is the geographic latitude of the telescope (37° 13' 55'' 7 in the case of OVRO), 5 is the 
declination of the object, and H is the hour angle. For objects transiting north of the zenith 
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(6 > A), ip = 180° at transit, while ip = 0° at transit for 5 < A. The position angles of the reference 
beams relative to the cluster are ±90° away from the parallactic angle at the given time. 

The reference beams are (very nearly) symmetric, so that we can fold the parallactic angle tp 
into the range (—90°, +90°) by adding or subtracting 180° with no loss of information. We define 
the principal parallactic angle (ip p ) in this way. It is desirable for the observing time on a cluster 
to be as evenly distributed as possible over a large range of principal parallactic angle, since this 
makes it easier to identify contaminating sources in the reference beams. However, clusters which 
transit near the zenith spend the majority of their time at the extremes of their ip p range. The 
7'35 beam FWHM at the switching separation of 22fl6 subtends Aip p ~ 19°; this is the effective 
'resolution' in parallactic angle. 

3.2. Lead- Trail Referencing 

The first and second levels of switching produce a FLUX measurement which is a "double- 
difference": the difference between the signal in a central 7'.35 FWHM beam and the average of 
two reference beams 22fl6 on either side of the central beam. However, even with this differencing, 
offsets in the FLUX levels remain at the ~ 100 level. These offsets change as the telescope 
tracks a cluster, therefore we must impose another level of differencing to remove this systematic 
effect. 

The third level of switching involves observations of LEAD and TRAIL control fields far from 
the cluster centers. The LEAD and TRAIL fields are offset by ±5a in right ascension from the 
cluster. Observations of the LEAD, MAIN, and TRAIL fields are separated by 5a m in time so that 
the telescope tracks through the same azimuth and zenith angles for each field. This referencing 
was used in the Coma observations of Herbig et al. (1995). From each MAIN FLUX, we subtract 
the average of the LEAD and TRAIL FLUXes to form the referenced field M - [L + T]/2 . This 
level of switching corrects for offsets dependent upon telescope orientation such as curvature in the 
ground spillover. We generally choose 10 m < 5a < 25 m , which is large enough so that the efficiency 
of the observations is not greatly reduced by frequent slewing, but small enough that temporal 
variations in the offsets are not severe. Though essential for the removal of systematic errors, the 
three-level switching technique employed greatly reduces the efficiency of our observations. Only 
one-sixth of the observing time is actually spent on the clusters, while the remainder is spent on 
the reference arcs and control fields. 

The criterion for selecting 5a is that both LEAD and TRAIL fields should be free of confusing 
radio sources bright enough to adversely affect the observations. We avoided sources found in the 
4.85 GHz Green Bank survey (Gregory & Condon 1991). For clusters with 5a = 15 m a typical 
scan consisted of six FLUXes on each of the LEAD, MAIN, and TRAIL fields. Each FLUX had a 
total integration time 4rj = 92 s . Fewer FLUXes could be completed per 15 m scan near the zenith 
because of the increased slew time there. The cluster A478 was observed with 5a = 19 m , which 



- 9 - 



allowed ten FLUXes per scan with integration time 4rj = 68 s . Approximately three minutes per 
scan were allocated for calibration and slewing. To ensure complete parallactic angle coverage 
for all three fields, observations were scheduled in different LST blocks. These different schedules 
filled in gaps that would have been created if the telescope had been slewing at the same LST in 
every schedule. The pointing positions (J2000) used for the cluster LEAD, MAIN, and TRAIL 
fields in the observations are listed in Table 

3.3. Calibration 

Calibration of the telescope and receiver was conducted regularly during the period from 
1993 to 1995, with the most exhaustive calibration measurements taken in 1995. We derive our 
calibration scale from the 1995 data and apply it to our SZE data taken in 1993 - 1994 using 
observations of standard sources. 

The telescope records the differenced power measurements for each FLUX in units of counts. 
In addition, the power from a calibrated noise diode is measured once at the end of each scan to 
normalize the changing gains in the system. The noise diode is calibrated in antenna temperature 
units by comparing its power output with hot (300 K) and cold (77 K) loads filling the beams, 
and in flux density units (1 Jansky = 10 -23 erg cm -2 s _1 Hz -1 ) by comparing against 'standard' 
radio sources such as Jupiter, Mars, and DR21. 

In addition to the waveguide switch that switches between the ANT and REF horns, the 
5.5-m receiver has an additional switch in each of the ANT and REF arms that allows us to switch 
between the horn and a cold load as inputs (SKY and LOAD). The value of the CAL is measured 
with both arms looking at the loads, since this improves the stability of the measurement. 

For each FLUX measurement, we assign a CAL value by linearly interpolating between the 
CAL measurements made before and after each FLUX. However, we do not use a CAL value 
if its standard deviation is greater than 5% of its value. This criterion rejects only about 1% 
of the CALs. We also do not use data if the difference between the bracketing CALs is greater 
than 5% of the mean value. Each FLUX and SD (standard deviation) measurement is divided 
by its assigned CAL value and multiplied by the appropriate temperature equivalent T ca [ to yield 
antenna temperature. 

The receiver system is slightly nonlinear. This nonlinearity has been measured and shown 
to be stable. Herbig (1993) has modeled the nonlinearity by assuming that the observed output 
power difference and the true input power difference are related by a linear function of the 
observed power. Because the CAL measurements are made with the load switches in the LOAD 
position, while the FLUX measurements are made with the switches in the SKY position, a 
correction must be made for the nonlinearity difference between the two power levels. We use the 
average total power, which is recorded once during each scan with the antenna looking first at the 
sky and then at the load, to correct the CAL value to the power level of the FLUX measurement. 
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Since the power levels change somewhat over the course of the observations, the correction must 
be applied to each FLUX individually. For the dataset as a whole, over a range of zenith angles 
and atmospheric noise levels, the correction factor ranged from 0.96 to 1.00. 

The aperture and beam efficiencies of the antenna are determined by comparing the flux 
densities of 'standard' sources to the measured antenna temperatures. Measurements of Jupiter, 
Mars, DR21, NGC7027 and 3C286 were made with the 5.5-m 32 GHz system during the period 
from February to June 1995. The flux density scale adopted for the 32 GHz observations are 
based upon a physical temperature of Tj = 144 ±8 K (Wrixson et al. 1971) for Jupiter at 
this frequency. This value encompasses the various measured temperatures for Jupiter, which 
range from 137 K to 153 K. DR21 and Mars are used as secondary calibrators. DR21 is an 
H II region in the galactic plane, and the 5.5-m measurements are confused by emission in 
the reference beams. The measured flux of DR21, relative to Jupiter assuming Tj = 144 K, is 
Sdr21 = 19.4 ± 0.3 Jy, in agreement with 19.4 Jy measured by Aller (private communication), 
though not with the 18.2 Jy value given by the Baars scale (Baars et al. 1977). We also measure a 
temperature of Mars relative to Jupiter giving Tm = 180 ± 3 K. The flux densities of the planetary 
nebula NGC7027 and of the radio galaxy 3C286 were also measured relative to DR21. We find 
S3C286 = 1-92 ± 0.06 Jy for S DR2 i = 19.4 Jy. Note that NGC7027 is expanding, and during Feb 
- June 1995 SjV7027 = 5.10 ± 0.14 Jy. The scale is listed in Table ||, both relative to Jupiter and 
using the adopted absolute value Tj = 144 K. 

The equivalent flux density and thermodynamic load temperature for the CAL diode are 
related by 

Sc^ = J±_ = 116.8 Jy/K 
T ca i A p rj A rj A 

where the physical area A p of the 5.5-m telescope is 23.64 m 2 . Using the Tj = 144 K flux density 
scale and measurements of the standard sources and hot and cold loads in April 1993, June 1993, 
October 1994, and April 1995, we calculate the aperture efficiency of the 5.5-m telescope to be 
r/A = 0.497 ± 0.007 (4.255 ± 0.014mK Jy" 1 ). The CAL temperature T cal corresponds to the 
equivalent Rayleigh-Jeans temperature increase in a blackbody filling the entire telescope beam. 
Because the power output of the CAL diode has been found on occasion to vary on timescales of 
months, the equivalent temperature or flux density of the CAL must be determined for the time 
period of calibration. During the SZE observing sessions from July 1993 to April 1994 the ratio 
of DR21 to the CAL was measured to be 0.1542 ± 0.0022. The stated uncertainty is the standard 
deviation about the mean for 12 observations spread over this interval and should reflect the range 
of possible variation in the CAL over this period. From these measurements we deduce that the 
CAL was stable to at least 1.4% rms. Assuming a flux density of 19.4 Jy for DR21, we adopt a 
constant value of S ca i = 125.8 ± 1.8 Jy, corresponding to T ca i = 535 ± 11 mK. After normalizing 
our FLUX values by the CAL, we multiply by T ca i in order to compare the measurements with 
the expected thermal noise. 

The relationship between the flux density and measured antenna temperature of a source 
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filling the main beam of the telescope is 

fr = fr On& = 161±4Jy/K (8) 

for the the measured main beam solid angle 0, m b = (5.12±0.14) x 10 -6 sr and observing wavelength 
A = 0.938 cm. Note that (0) and (§) differ by 

— = = 1-38 ± 0.04. (9) 

where r]s is the fraction of the entire telescope beam in the forward main beam. From the 
aperture efficiency and main beam solid angle we derive a beam efficiency ijb = 0.684 ± 0.021. 
A radio source of flux density S ca i is equivalent to a source filling the main beam with uniform 
Rayleigh-Jeans brightness temperature 

T* ca i = S cal = — = 782 ± 24 mK. (10) 

This is the appropriate CAL value to use for calibrating the astronomical signal in the main beam. 
The FLUX measurements, already scaled by T ca i, must be divided by the beam efficiency t\b to 
effectively scale by T* al . 

The FLUX and SD values must also be scaled by a factor which accounts for atmospheric 
attenuation. Atmospheric attenuation reduces the intensity by a factor ~ e -asece z ^ w h ere Q z { s 
the zenith angle and a is the optical depth. In good weather at v = 32 GHz, a ~ 0.04 (assuming 
a ~ Tatm/Tphyg ~ 11K/270K). The estimated rms variation in the atmospheric optical depth is 
°a 55 0.02, based upon water- vapor radiometry data obtained at OVRO. 

In summary, the raw FLUX and SD values are first divided by the CAL value, then scaled by 
the CAL temperature T ca i. After editing, the measurements are scaled by the correction factor : 

a sec Bz 

k = = 1.46e a04sec ^. (11) 

VB 

At this stage, the measurements are in degrees Kelvin, and are equivalent to the differences in 
temperature between two high-temperature uniform blackbody emitters (calculated using the 
Rayleigh-Jeans formula) filling the main beams of the ANT and REF horns that would produce 
the observed power differences. 

The total calibration uncertainty is the quadrature sum of the uncertainties in: (1) the 
temperature of Jupiter Tj (5.6%), (2) the ratio of our measured flux of DR21 relative to Jupiter 
(1.6%), (3) the ratio of the CAL S ca [ to DR21 during the course of the observations (1.4%), (4) 
the main beam solid angle converting flux density to brightness temperature (2.7%), and (5) the 
atmospheric attenuation variation (2.0%), giving a total calibration error budget of 6.9%. The 
dominant uncertainty is the absolute brightness temperature of Jupiter. The relative flux density 
scale, (2) and (3), is accurate to 2.1%, while the temperature scale, (2), (3) and (4), is accurate to 
3.4% (without atmospheric correction). For improvement in the SZE measurements beyond what 
is reported here, a more accurate absolute flux density scale at 32 GHz is required. 
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3.4. Data Editing 

Much of the data taken during the afternoon or when the weather was bad was clearly 
unusable, and a method of editing the data without introducing any systematic biases had to be 
devised. The two general editing methods used were editing based on the SD values and filtering 
based on the standard deviation of the switched difference measurements. 

Before applying these methods, we removed FLUXes for which the elapsed time exceeded the 
expected duration by more than 4 seconds. Excessive durations could be caused by extremely 
high winds or tracking problems which would corrupt the data. In addition, on several occasions 
problems arose with the zenith angle encoder which caused the drive to fail, and resulted in a few 
extremely long FLUXes during which the telescope was not tracking the source. 

The procedure for thermal editing is based upon the expected thermal noise in a FLUX 
measurement. The thermal standard deviation is given by: 

V4Tj Au \jTi lAv 

where 4r« is the total integration time of the FLUX, Av = 5.7 GHz is the bandwidth of TRW 
amplifier, and the numerical factor arises because of the double switching inherent in the FLUX 
measurement and the FLUX definition of (||). The system temperature T sys is given by 

T sys = T + T atm sec 9 Z - (13) 

T a tm is the atmospheric emission temperature per airmass. Included in To are contributions from 
the receiver, the ground, and the 2.726 K cosmic background (Mather et al. 1994), 

T = T rx + T gnd + 2.726 K (14) 
T rx = 33.0 ± 1.6 K 

Tgnd ~ 8 K 

T atm ~ 10 K. 

At one airmass for r, = 20 s , this gives Oth = 0.16 mK. Even with the overhead from the three-level 
switching, we would expect to integrate down to a noise level of around 36 /iK in the referenced 
field M — [L + T]/2 in about 3 hours of total integration time. In practice, additional atmospheric 
noise not removed by the double differencing increases the noise significantly above this level. 

We first wish to reject blocks of data which are clearly contaminated by bad weather. Our 
method is a modification of the method developed by Brandt (1992), in which one searches for 
blocks of good data which contain a specified number of points, typically 10 to 25, within a limited 
time range, typically 2 hours. Sliding buffers are moved over the dataset, and any point that is 
contained within a "good" block of data is accepted. 

For each point i in the dataset, we construct the "test" statistic 

ti = Xi/a t h,i- (15) 
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We have divided the value by the corresponding thermal noise level ath,i to account for the 
increased scatter expected at high zenith angle, where the telescope is looking through a longer 
column of air, and to allow comparison of data points with differing integration times. The value 
X used as the statistic can be FLUX, SD, SW1 or SW2. 

The test values are placed into "buffers" of N consecutive points, for which the mean and the 
standard deviation about the mean are calculated. For buffer j, 

1 j+N-l 

tj = - Yl u ( 16a ) 

i=j 



i j+N-l 



i=j 



1/2 



(16b) 



Each test buffer contains N points, restricted to be taken over a maximum span of not more than 
2 hours. 

The filter is effected by examining all buffers that contain a given data point i. A data point 
i is rejected if there exists no buffer j of the chosen length N containing the given point for which 
the standard deviation Sj is less than a chosen cutoff value s max and for which the mean tj is less 
than a limit t max . We used either the mean or the standard deviation of the test values X in the 
filter, but not both at once. We designate these filters as "meanXfiV, t ma a;]" or "sigXfiV, s ma x]" ■ 
In addition, W6 Cciii simply throw out discrepant values with t% !> tmaxi 

with filter designation 

ll X[t max ]" . This brute-force rejection is only useful for the SD, where it can remove single points 
with large error bars that escape previous filters. In these filters, the FLUX values themselves are 
not used to avoid bias of the results. 

Through experimentation upon the data, we have developed a filtering sequence that produces 
good results and is robust to slight changes in the filter parameters: 

1. Filter on mean of SD, with N = 25 and t max = 2.0 (meanSD[25, 2.0]). 

2. Reject points with high SD, using t max = 2.5 (SD[2.5]). 

3. Filter on standard deviation of SW1, with N = 10 and s max = 7.5 (sigSWl[10, 7.5]). 



After our filtering, we are left with a "clean" distribution of FLUXes, although a few FLUXes 
with discrepant values from the mean may remain. These points can escape our culling procedure 
if they are due to stationary structures in the far-field atmosphere, or, of course, spurious noise or 
interference that happens to mimic the switching scheme of the FLUX (unlikely, but possible). 
The removal of "bad" FLUXes which show no ill effects in the respective SW1, SW2, or SD 
without biasing the data is a tricky business. A conservative approach would be to reject all 
FLUXes that are part of a scan on a given field for which the standard deviation of the FLUXes 
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about the scan mean is above some limit — this should be equivalent to a noise- level edit. In 
practice, we will use the scan standard deviations to down-weight this data during the analysis 
(see next section). 

Our approach, one we have adopted with success in the past (Readhead et al. 1989, Myers et 
al. 1993), is to throw out the FLUXes on the tail of the FLUX distribution for a given field using 
an iterative procedure to ensure that the resulting post-edit distribution is not skewed by applying 
too stringent a cutoff. This rejection procedure we designate as "rejX[a]". For the cluster data, we 
reject iteratively the referenced temperatures ATmlt (see below) with a cutoff of 4cr (rejMLT[4]). 
Note that for reasonable cutoffs (> 3 a), this procedure will not introduce a bias into the data, 
since outliers are rejected one at a time and after each iteration the mean and standard deviation 
are re-calculated. 

Discussion of the effects of different editing and filtering parameters will be presented in the 
next section. 



3.5. Analysis of SZE Measurements 

Software was designed to perform the subtraction of the reference LEAD and TRAIL fields 
from the MAIN field. The program first multiplies the filtered and edited data by the correction 
constant k from equation (fij 



ATj = kFLUXj a i = KSD i (17) 

to convert the power differences into Rayleigh-Jeans temperature differences in degrees Kelvin 
between the ANT and REF main beams. The separate FLUXes from the adjacent MAIN (M), 
LEAD (L), and TRAIL (T) scans are then matched based on their proximity in the azimuth 
and zenith coordinates. The FLUXes in the MAIN fields are matched with the closest FLUXes 
in the adjacent LEAD and TRAIL fields that have not already been matched to other MAIN 
points. Matching LEAD and TRAIL measurements are required for each MAIN — if no reference 
FLUXes are found in the LEAD or TRAIL that are closer than 5' (a substantial fraction of the 
beamwidth) from the MAIN field, the MAIN FLUX is discarded. 

For each triplet, we calculate the "referenced" differences 

ATmlt = AT M - ^{AT L + AT T ) (18a) 
ATm-l = ATm-AT l (18b) 
AT M _ T = AT M -AT T (18c) 

ATjTf = 1 -{AT L + AT T ) (18d) 
ATl-t = AT L -AT T , (18e) 
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Gaussian error propagation is used to compute the statistical uncertainties in these differences. 
For example, 

^MLT = 4r + 4^i + 4)- (19) 

The referenced points are then weighted and the statistics are formed. For each of the N 
points with M - [L + T]/2 , M - L , M - T , etc., we form AT ± e where 

I N 

AT = ^X>; AT i ( 20a ) 

, 2 = (20b) 



a 2 



N 1 N 



N — 1 \\2 ■_ i 



iV 

Wi = Y. w i ( 20d ) 

j'=i 

N 

w 2 = J2 w l ( 20e ) 

3=1 



The weights are formed from a combination of the individual errors aj (e.g. equation |l7|) and 
the standard deviation of the points within scans cr sc , 

wj 1 = a + a\ a 2 + a 2 a 2 c j. (21) 

The standard deviation SD recorded with each FLUX generally underestimates the actual error 
by a factor of 2 or 3. This is a result of atmospheric fluctuations on a timescale longer than 
the duration of the FLUXes. A better estimate of the true error is the scatter of the FLUX 
measurements within a scan. An appropriate weighting is a\ = a<i = 1 and ao = 0. The inclusion 
of the individual SDs a\ ^ guards against anomalously low scan standard deviations causing 
very large weights, while the weight will normally be dominated by the larger scan standard 
deviations cr sc j. 

Fig |l] shows our results for A478, A2142 and A2256 using twenty different filtering methods. 
The referenced measurements ATmlt are plotted. Both unweighted (u: ao = 1, ai = a-i = 0) and 
weighted (w: a± = a% = 1, ao = 0) results are shown. The methods are enumerated by filtering 
and outlier editing schemes. There were five filtering methods chosen: 



Method 1: No filtering 

Method 5: SD[3.0] 

Method 9: meanSD[25, 2.5], SD[3.0] 



-16- 



Method 13: meanSD[25, 2.0], SD[2.5], sigSWl[10, 7.5] 
Method 17: meanSD[25, 1.5], SD[2.0], sigSWl[10, 6.0] 

In addition, for each filter, four different outlier rejection schemes were tested: 

Method n+0: No outlier editing 
Method n+1: rejMLT[5] (5a rejection) 
Method n+2: rejMLT[4] (4a rejection) 
Method n+3: rejMLT[3] (3a rejection) 

These filter methods 1-20 are roughly in increasing order of fraction of the data rejected. A 
total of 802.5 hours of MAIN, LEAD and TRAIL data were passed with no filtering, editing, or 
weighting (Method lu), while only 351.7 hours of data were accepted for Method 20u. 



For the weighted means, the form of ( |20a| ) suggests an effective number of points 

Neff - W - ^L^f (22) 

The effective fraction (~ N e ff/N) of accepted data ranges from 35% (Method lw) to 27% 
(Method 20w) with little variation between outlier editing modes for each filter method. Note the 
smaller variation between the weighted methods; the weighting largely takes care of the editing 
by down-weighting bad stretches of data. The unedited and unweighted data point (Method lu) 
is not visible in Fig |], as it is heavily corrupted by bad data and is off-scale. The corresponding 
weighted point (Method lw) is consistent with the other data points. 

We adopt method 15w: meanSD[25, 2.0], SD[2.5], sigSWl[10, 7.5] filtering, rejMLT[4] outlier 
editing, and weighting a\ = a% = 1, ao = (marked in Fig [[]). Because none of the results differ 
significantly with respect to the statistical error bars, we are confident that even fairly large 
differences in the filtering method have negligible effect upon our results. We conclude that our 
automatic editing procedure is reliable and robust. 



4. SZE Results 

In all three clusters, A478, A2142, and A2256, we find significant detections of a microwave 
decrement that we attribute to the SZE. Table || lists the measured values of AT for these clusters. 
In addition to the fully referenced difference M — [L + T]/2 , and the single-referenced differences 
M — L and M — T , we also list the averages for each of the MAIN, LEAD, and TRAIL fields 
separately, as well as the weighted difference and average of the LEAD and TRAIL. If there 
is no significant source contamination, we expect the averaged L — T to be consistent with 
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zero, provided that our switching technique is subtracting all significant ground spillover effects. 
However, if ground spillover or atmospheric emission varies with a timescale shorter than the 
time between the MAIN, LEAD, and TRAIL scans, the L — T may be inconsistent with zero, 
but this effect could be expected to average out over the course of an observing season. Source 
contamination in the control fields will generally result in non-zero values for L — T and will 
contaminate the average (L + T)/2 . 

The measurements of the SZE as a function of principal parallactic angle ip p are presented 
in Fig ^ - Fig ||. The upper panels of each plot show the individual MAIN, LEAD and TRAIL 
binned averages, while the lower panels show the referenced M — [L + T]/2 and L — T . The 
scatter of these referenced values with parallactic angle can give an estimate of possible reference 
arc contamination. 

For A478 we find AT = —375 db 24 /zK with 64 total hours of integration time, a 15.6cr 
detection. Fig || shows AT binned by parallactic angle ip p . The close agreement of the LEAD and 
TRAIL over the entire range of ip p suggests that there is no significant source contamination of the 
control fields. Although the M — [L + T]/2 is nearly constant with tp p , with x 2 = 1-70 about the 
mean in the 5 bins, the individual fields in the upper figure show variations of ~ 400 fiK resulting 
from ground spillover. This demonstrates the necessity of our LEAD and TRAIL referencing. The 
LEAD and TRAIL difference L — T is consistent with zero, and has x 2 = 6.87 about zero in the 
5 parallactic angle bins, which is significant only at the 79% level. Thus, we conclude that there is 
no evidence for source contamination in the A478 data. 

A microwave decrement of AT = —420 db 19 /iK is found in A2142, a detection significant 
at the 22ct level. There were 85 hours of usable integration time on this cluster. As shown in 
Fig |3|, the LEAD and TRAIL track each other with an average offset of 188 /iK. There is a mean 
L — T difference of —66 ± 20, which is marginally significant at the 3a level. In Fig |3| we see 
a feature in the L — T at parallactic angle ip p ~ —20°. This is consistent with a source in the 
LEAD reference arc centered at ifj p = —21° (see next section). The x 2 ls 3-14 about zero in the 
4 parallactic angle bins, although since most of the data is at the ends of the arcs, the effective 
number of degrees of freedom ( p2[ ) is only 2.7. The significance of this chi-squared value is 67%. 
The referenced M — [L + T]/2 has a x 2 °f 2-1 about the mean for 1.5 effective degrees of freedom, 
which is significant only at the 75% level. 

We measure AT = —218 ± 14 fj,K for A2256, a detection significant at 15.6o\ There was 310 
hours of integration time on A2256. The MAIN, LEAD, and TRAIL data are shown binned by 
parallactic angle in Fig |j. Because A2256 is circumpolar, the magnitude of the average LEAD 
and TRAIL is somewhat smaller for this cluster as it traverses a more restricted range in azimuth 
and zenith angle and is not as strongly affected by changes in ground spillover as are A2142 
and A478. The overall average difference between LEAD and TRAIL is again consistent with 
zero, —22 ± 15 /iK. However, when the data is binned in parallactic angle as in Fig ||, we find 
X 2 = 7.88 versus zero for 4.3 effective degrees of freedom, which is significant at the 88% level. 
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After referencing, the x 2 of M - [L + T]/2 data versus the mean is 5.17 for 3.3 effective degrees 
of freedom, or 81% significance. 

In all three cases, the MAIN field is clearly showing the SZE decrement relative to the LEAD 
and TRAIL, while all three track together with parallactic angle. The variation of the LEAD 
and TRAIL average offset with ip p shows the importance of the extra referencing to remove the 
differential ground spillover component. This referencing was also performed using a slightly 
different procedure in the Coma observations of Herbig et al. (1995). 

The dataset is also robust with respect to separating into different time periods and different 
times of day. No trends are seen that are significant compared to the statistical uncertainties. 

After referencing, we have rms uncertainty levels of 14 /^K-24 /j,K. On these angular scales 
at this fluctuation level, anisotropics in the cosmic microwave background radiation itself can be 
expected to be detectable. Intrinsic fluctuations on this scale are expected to be in the range 
5 x 1(T 6 £ AT/T £ 2 x 1(T 5 (14 ^uK-55 fjK) in the most popular models (eg. Bond et al. 1994). 
Our instrumental filtering can be expected to reduce the CMB anisotropy signal somewhat, as 
will the smearing with parallactic angle of the reference beams. The clearest indicator of CMB 
fluctuations would be significant L — T differences; the fact that we have no clear detection of a 
LEAD - TRAIL difference, except possibly in A2142, indicates that the true anisotropies are not 
much greater than the predicted range. In a separate observing program, we have conducted a 
microwave background anisotropy experiment using the 5.5-m telescope and the same instrumental 
configuration, the results of which will be reported in an upcoming paper. If the background 
fluctuations are indeed in the expected range, then pushing the SZE on these scales to much lower 
noise levels will not be possible using single-frequency measurements such as ours. Note, however, 
the SZE in more distant and hence much smaller angular-sized clusters will not be so badly 
affected, as the background fluctuations on smaller scales are expected to be considerably smaller. 

A significant contribution to systematic error in the SZE measurements is foreground 
contamination by galactic and extragalactic emission. At this frequency and angular scale we 
believe that galactic dust and free-free emission are not likely to be major contaminants. However, 
synchrotron emission by discrete extragalactic sources is known to be a significant problem. 

5. Source Contamination 

An unfortunate aspect of centimeter wavelength observations of the SZE is that they must 
contend with the presence of radio sources that can mimic or hide the effect. High resolution 
radio maps must be made at or near the frequency of observation to deal effectively with this 
problem. Unfortunately, many radio sources are also variable, and if a cluster is contaminated by 
such sources, be they field object or associated with the cluster, simultaneous observations must 
be made on different telescopes. This was not feasible for the present work; however, previously 
published observations of the fields of many of our clusters suggest that our results are not greatly 
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affected by source contamination. Note that in this respect interferometric SZE observations are 
superior — the longer baselines provide the simultaneous high-resolution information necessary 
for source identification, and if the baselines extend far enough compared to the shortest spacings, 
a clean subtraction of the interfering sources can be made. However, with single-dish data, we 
have not this luxury. 

Radio observations of adequate resolution for computing corrections to SZE measurements 
have been published for many of the clusters in our sample. However, these maps extend 
only about 20' from the cluster centers, so the possibility exists that there are unrecognized 
contaminating sources in the reference arcs. For the LEAD and TRAIL fields, we use data from 
the 1987 Green Bank survey at 4.85 GHz (Gregory & Condon 1991), which lists sources down 
to a flux limit of 25 mJy. Unfortunately, the declination range of this survey does not include 
A2256, and the radio environments of the LEAD and TRAIL field for this cluster are at present 
unknown. However, the excellent match between the LEAD and TRAIL fluxes and the stability of 
the referenced ATmlt with ip p give us confidence that contamination is not a serious problem in 
A2256. Many of the clusters in our sample fall below the southern declination limit of the Green 
Bank survey, and for these we obtained source positions from preliminary results of the southern 
Parkes-MIT-NRAO (PMN) survey (Griffith & Wright 1994). 



5.1. 40-m Observations at OVRO 

The published observations of discrete radio sources have been carried out at frequencies 
considerably lower than our observing frequency of 32 GHz. For this reason, we observed the 
contaminating sources from the Green Bank and PMN surveys with the OVRO 40-m telescope 
at 18.5 GHz. We observed sources within 30 m in right ascension and 32' in declination from the 
cluster center, in order to cover the possible locations for the LEAD and TRAIL fields. Sources 
near the cluster centers were also identified from higher resolution radio maps. Images of A478 
and A2142 at 2.7, 4.75, and 10.7 GHz were available (Andernach et al. 1986), as well as for A2256 
at 610 and 1415 MHz (Bridle et al. 1979). 

At a frequency of 18.5 GHz, the 40-m telescope has a beamwidth of 2' FWHM. Because 
of its larger collecting area and much higher sensitivity to point sources, the 40-m can measure 
any sources that would be bright enough to affect 5.5-m SZE observations in a relatively short 
integration time. Observations of these confusing sources were carried out during the period 
between 20 November 1993 and 24 January 1994. Many observations were repeated several weeks 
after the initial observations in order to gauge variability. 

The desired sensitivity level was achieved in around 400 s of integration time on each source. 
On the 40-m telescope, we use the same double differencing procedure in measuring FLUXes as on 



the 5.5-m (§ 3T). Pointing was checked before each scan on a source. Calibration was performed 



by observing DR21, which has a flux of 19.2 ± 0.7 mJy at 18.5 GHz. Because of the large size of 
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the 40-m telescope, physical deformation of the dish causes the gain to vary with elevation. Long 
tracks on 3C84 were used to determine the zenith-angle dependent gain corrections. We estimate 
the error for calibration and gain corrections at ~ 6%, similar to those for the 5.5-m telescope. 

The 40-m 18.5 GHz measurements are listed in Table S, along with the lower frequency flux 
density measurements (Si ow ) at 4.85 GHz, 2.7 GHz and 10.7 GHz obtained from the literature. 
Where the 18.5 GHz measurement was not significant at the 3a level, the 3a upper limit on the 
flux density is listed. The deduced spectral index a for each source, between the lower frequency 
and 18.5 GHz, are given, where a power-law spectrum is assumed 

S cx v a . (23) 

The extrapolated flux densities S32 are listed in the last column. These were calculated using the 
values for a listed in the table, or from upper limits where appropriate. 

Because the 18.5 GHz frequency of our source measurements is significantly lower than the 32 
GHz observing frequency of our SZE data, and we have no bracketing measurements at higher 
frequencies, one should conservatively assume the 18.5 GHz measurements themselves as an upper 
limit to the 32 GHz source flux densities. It is possible that the 18.5 GHz emission is dominated 
by flat-spectrum compact components in the radio sources. From the numbers in Table || we 
see that this leaves us with an overall factor of two uncertainty in the source corrections to be 
applied where a detection at 18.5 GHz was made. Only better measurements of these sources at 
frequencies bracketing 32 GHz will allow accurate source subtraction to be made. 



5.2. Corrections to SZE Measurements 

The radial distances and position angles (relative to North through East) of these 
contaminating sources relative to the field centers are given in Table |. If the source lies within 
the reference beam arc 22fl6 away from the central beam, then the corresponding distance from 
the closest arc center r and the parallactic angle vfj p of the closest approach are also listed. The 
implied equivalent SZE temperatures AT a are calculated using 6.22 mKJy -1 , since the SZE 
measurements have already been converted to main beam temperatures (see equation |||). There 
is a 6.6% overall conversion uncertainty that includes all of the flux density scale uncertainties 



(see § 3.3), although this is less than the uncertainty in the source flux density extrapolation to 
32 GHz, and thus has not been applied to the numbers in Table |6|. Finally, correction for the 7'.35 
Gaussian beam pattern produces the expected maximum contributions AT max printed in the last 
column of the table. 

In the MAIN fields, if a source falls within the central part of the main 5.5-m beam, it will 
tend to cancel out the SZE decrement. Sources in the reference arcs are subtracted and thus 
mimic a decrement, but they will only be observed at certain parallactic angles. Such sources are 
recognizable by their signature on a plot of the SZE decrement vs. parallactic angle ip„. Sources 
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in the beams of the LEAD and TRAIL fields will have the opposite effect compared to those in 
the MAIN field, with levels reduced by a factor of two due to the averaging ( M — [L + T]/2 ). 

Because each FLUX measurement is taken at a given parallactic angle, and enters into the 
final mean with its individual weight, it is necessary to subtract the effects of sources in the 
reference arcs on a point-by-point basis. In practice the LEAD, MAIN and TRAIL measurements 
for each referenced measurement are adjusted by the values determined using the positions and 
AT a from Table |(| Only those sources with significant 18.5 GHz detections were used. In Table 
we list the unweighted mean corrections AT^ applied to the data points in the contaminated 
fields. The center beam corrections AT on and maximum reference arc corrections AT re y (for the 
parallactic angle where the source is closest to the center of the arc) are also given. 

The corrected SZE results are shown in Table ||. The second column lists the source 
contributions AT5.5 m (src) computed using the actual data weighting in the cases of A2142 and 
A2256, and the 3<r limit for A478 (source A478.1). The corrected values AT5.5 m (corr) are listed 
in the final column, with the correction values themselves added in quadrature as an uncertainty. 
For A478, the la limit on the contribution of source A478.1 was used as the uncertainty. The 
measurement of the SZE in the weakest cluster A2256 is most adversely affected by the correction 
uncertainties, with the statistical standard error increasing from 6.4% to 12%. The error bars on 
the A478 and A2142 are not as strongly affected due to the larger relative decrements. However, 
in all three cases the corrections applied were similar in magnitude to or larger than the purely 
statistical measurement uncertainties, and thus contribute significantly to the error budgets. More 
accurate SZE measurements will require better source measurements, necessarily contemporaneous 
with the SZE observations to deal with possible variability in the foreground sources. 



6. The SZE, Baryonic Mass, and H 

The SZE is proportional to the Compton y-parameter 

oo 

kT f 



y 



f _^ OT n edC> . (24) 
J m e c z 



For convenience, we will use cylindrical coordinates (R, 4>, £) (R, £ in Mpc) centered upon the 
cluster with £ along the line-of-sight. Then 



oo 
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y(R,0)= [ -^n e (R,<l>,OT e (R,<f>,C)dC. (25) 
J m e c z 
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For the small angles considered here, R = D a 9 for angular diameter distance D a to the cluster. 
For g = 1/2 and H = 100/ikms" 1 Mpc -1 , which we will assume throughout this papern, 



D a = 6000 (1 h- 1 Mpc. (26) 

(i + zy 

The on-sky intensity differences are measured in units of antenna temperature AT a , which is 
the equivalent temperature difference in the Rayleigh-Jeans limit 

AI„ = ^AT a . (27) 

Using the standard formulae for the fractional change in the intensity of the thermal background 
in the non-relativistic limit (Sunyaev & Zeldovich 1980), we get the frequency dependence of the 
measured change in antenna temperature of the microwave background due to the SZE: 

AT a x 2 e x ( , x \ 

— i = y - _ x coth - - 4 , (28) 

T cmb y (e*-l) 2 \ 2 J> K J 

where x is the dimensionless frequency 

hu v , 

x = = (29) 

kT cmb 56.80 GHz' v ' 

We use the COBE FIRAS value for the microwave background temperature T cmb = 2.726 ±0.010 K 
(Mather et al. 1994). At the 5.5-m observing frequency v = 32 GHz (x = 0.563) this is 

AT 

_«= -1.8971/. (30) 

-L crab 

There are several possible corrections to the expressions (|28|) and hence (|30|). In addition to 
the thermal SZE, there is a kinematic effect due to the peculiar velocity of the cluster (Sunyaev and 
Zeldovich 1980). At 32 GHz, a peculiar velocity of 300 km s -1 for a cluster with kT e = 7.5 keV will 
produce a change in the SZE intensity of only 2%, and can safely be ignored in these calculations. 

Another factor not accounted for in our expression for y is the relativistic correction to (pq), 
which was derived in the non-relativistic limit. Rephaeli (1995) has calculated the corrections 
for the low optical depths (r < 10~ 2 ) and mildly relativistic electron temperatures (kT e ~ 5 - 
lOkeV) appropriate to our clusters. For these parameters, Rephaeli finds corrections of around 
+3% ± 0.3% at our observing x = 0.563, i.e., the SZE decrement is less pronounced in magnitude 
than what eq. (|28|) would predict. Because this is a systematic underestimation of y given an 
observed AT a , we use the corrected relation 

Xrel AT a /Q1\ 

Vmeas ~ ~ 1.897 T cmb 1 j 



3 At these low redshifts, the effect of the cosmology is purely kinematic and thus depends only upon go- To first 
order, AD a /D a ~ Ago z/2, which at our redshift limit of z — 0.1 makes a ±2.5% change in D a , and thus the derived 
h, for Ago = ±1/2. 
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where the relativistic correction factors Xrei for each cluster are given in Table 10, along with the 
Compton parameters y me as using the source-corrected AT a from Table ||. 

Herbig et al. (1995) carried out SZE observations and analysis of the Coma cluster using 
the same 5.5-m setup and calibration scale as was employed in our observations. With source 
corrections and calibration uncertainty included, they found a non-relativistic y-parameter of 
1/meas = (5.96 ± 0.99) x 10~ 5 . If we apply the relativistic correction 1.029 to the Herbig et 
al. measurement, and remove the 6.9% calibration uncertainty, we get y me as = (6.13 ±0.93) x 10 -5 . 



This is the value we have given in Table 1C. 



The observed SZE decrement is the true decrement modified by the telescope primary beam 
and the beam switching. The 5.5-m single beam pattern is well approximated by a circular 
Gaussian 

= d«| exp (32) 

with beam width 6 g = 3^12 ± 0f04 (7f35 ± OflO FWHM). The beam is less than 1.4% elliptical. 
Because clusters track through a range of parallactic angles, the slightly elliptical beam is rotated 
on-sky and therefore the effective beam is well represented by the geometric mean 9 g . 

The average Compton y-parameter in the Gaussian beam, on a line of sight offset at cylindrical 
radius R from the center at position angle </>, for small angles is given by 

2n oo 

y g (R, (f>) = j d6 j ' r dr^—^e- d2/2L sy(r, </> - 0) d 2 = R 2 + r 2 - 2Rr cos(0). (33) 



o o 



9 



Here L g = D a 6 g . The beam-switching can be evaluated as 



y sw (R,<t>;ip p ) = y g {RA) - ^y g {R-,<P-) - ^y g {R+,(f>+) (34) 

with 

R 2 ± = D 2 sm 2 (V?p + <p) + [R±D cos(V> P + <p)f (35a) 
J Rshi(/>T£>sinV>p 

tan</> ± = — — f-. (35b) 

R cos (p ± L> cos ipp 

The 5.5-m beam separation is 6 a = 22fl6, so D = D a 0£>, and ip p is the principal parallactic angle 
for the observation. For convenience, the position angle 4> is measured starting from the East so it 
is in the same orientation as the parallactic angle tp. 

We will most often use the integrals through the cluster center in a cylindrically symmetric 
model. The angular dependences are dropped, so y g (R,(f>) = y g (R) is now a function of angular 
radius from the cluster center only. The expressions (^3[) and (|34]) are abbreviated as 

oo 

y g = y g (0) = 2vr J r dr^e^ 2 ' 2L ly{r). (36) 
o 9 



and 



y sw = y sw (0) = y g (0) - y g {D). 



(37) 



6.1. Density models for the ICM 

The distribution of the X-ray emitting gas in galaxy clusters has frequently been modeled by 
an isothermal /3-model (Cavaliere & Fusco-Femiano 1976), also known as a modified isothermal 
King model. In this spherically symmetric model, the electron gas density n e is given as a function 
of the spherical radius r from the center of the cluster by 

n e = n {l + ^j , (38) 

where r c is the core radius of the cluster and no is the density at r = 0. For the present, we will 
consider mass distributions with circular symmetry in the plane of the sky, although it is easy to 
generalize to ellipsoidal profiles. 

Table |9| lists published parameters of the intracluster gas from various X-ray observations 
of A478, A2142, A2256, and Coma. The core radius 9 core determined from the X-ray surface 
brightness profile is listed, and has been converted into a linear size r c in h- 1 Mpc using (EBT). The 
overall temperature of the X-ray emitting gas T e is given as kT e in keV (10 s K = 8.61 keV). As 
determined from the X-ray emission, the central densities no are in units of /i 1//2 cm~ 3 . A value 
for (3 is given only if it has been determined by fitting the surface brightness profile. All errors are 
given as la, converted from 90% confidence (~ 1.645<r) in the literature if necessary. 

We now discuss the models for each cluster in detail. Because the model uncertainties are in 
most cases the dominant source of systematic error in the determinations of the baryonic masses 
and Hubble constant for this sample, we plan to make our own detailed analysis of the ROSAT 
data for these clusters to obtain more accurate models, and to understand the limitations of our 
particular method better. For now, we adopt the models presented in the literature and proceed 
with our analysis. 



6.1.1. A478 

Abell 478 has been shown to contain one of the largest cluster cooling flows with more 
than 5 x 10 11 h~ l M of X-ray absorbing matter within the inner 150 h~ l kpc and a total mass 
deposition rate of ~ 500 h^ 1 M Q yr _1 (Johnstone et al. 1992, Allen et al. 1993). The combined 
Ginga and ROSAT data measure the temperature kT e = 6.56 db 0.09 keV for the cluster isothermal 
component, and a temperature of kT e ~ 3keV within the inner 75 h~ 1 kpc. The ROSAT PSPC 
image of A478 (Allen et al. 1993) shows an axial ratio of ~ 0.8 to the inner (8 < 2'A) isophotal 
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contours. This would suggest that A478 is probably even more ellipsoidal than this, and should 
be kept in mind in the proceeding analysis. 

The presence of such a large cooling component to the cluster core medium makes modeling 
of the SZE from the X-ray emission difficult. Edge & Stewart (1991a) found T e = 6.8 keV, 
r c = 0.10 h -1 Mpc and n = 25.2 ± 2.8 x 1CT 3 h 1 / 2 cm~ 3 from EXOSAT observations. Allen 
et al. (1993) fit a central electron density to the isothermal (non-cooling) component of 
no ~ 9.55 x 10~ 3 /i 1 / 2 cm -3 assuming a core radius of r c = 0.125 Mpc and a King profile 
09 = 2/3). 

We have adopted the Allen et al. (1993) values for no, r c and ft The uncertainties on these 
quantities were determined empirically by comparison with our own preliminary analysis of the 
ROSAT data. This model is the one listed in Table |9[ 

6.1.2. A2142 

A2142 is the most distant cluster [z = 0.0899) and has the largest 2-10 keV luminosity in our 
sample. A2142 is also the second most luminous cluster in the Edge sample as a whole. Edge & 
Stewart (1991a) list A2142 as a cooling core cluster, and Edge, Stewart & Fabian (1992) derive a 
mass flow rate of 50-150 h~ l M Q yr _1 . 

In the compilation of cluster temperatures by David et al. (1993), A2142 is listed as 
having a Ginga temperature of kT e = 8.68 ±0.12keV. Abramopolous & Ku (1983) derive 
r c = 0.26 ± 0.01 h~ l Mpc and n = 6.97 ± 0.41 x 10~ 3 h 1 / 2 cm , where we have estimated the 

1/2 —3/2 

uncertainty in no through the relation no oc L x r c . The parameter (3 is fixed at unity in this 
model. For want of a better determination, we adopt /3 = 1 ± 0.3. 

6.1.3. A 225 6 

Although A2256 does not appear to have a central cooling flow (Edge, Stewart & Fabian 
1992), it does show significant substructure in the X-ray emitting gas. Briel et al. (1991) 
found evidence for a "merger event" in the ROSAT PSPC image of the cluster. Two surface 
brightness peaks were found in the cluster center with a separation of 3f5 (160 h^ 1 kpc), with 
some indication of differing temperatures. David et al. (1993) list an overall Ginga temperature 
kT e = 7.51 ± 0.11 keV, while the ROSAT PSPC data gives kT e ~ 2.0 keV for the cooler (NW) 
subcluster. Fits of the azimuthally averaged data excluding the secondary maximum to a modified 
isothermal King profile yielded C)1 = 4'.83 ± Of 17 and ft = 0.756 ± 0.013, while a fit to the 
secondary after subtraction of the primary smooth profile gave # Cj 2 = 4f3 ± 0f4 and ft = 1-1 ± 0.1 
and a peak surface brightness 82% of that of the primary. Briel et al. (1991) also analyze the 
distribution of 87 galaxies to find radial velocity dispersions of 1270 ± 127 km s -1 for the main 
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cluster region and 250 ± 123kms _1 in the NW subcluster. A relative systemic velocity difference 
of —2150 ± 259 km s -1 is found between the NW region and the main cluster — this difference is 
consistent with the infall velocity at 1 Mpc from a 10 15 M cluster. The pointing center used in 
our observations is approximately in the center, between the two components of A2256. 

Davis & Mushotzky (1993) examined Einstein IPC data and also found evidence for the 
merger, and derive 6 C = 6'0~^q® and [3 = 0.72^'^. By fitting elliptical isophotes, we find axial 
ratios of 0.6 near the center to 0.25 at a radius of 10' . Spectroscopy from the BBXRT (Miyaji et 
al. 1993) indicates a temperature of kT e = 4.6^'^ keV for the NW component. 

A more detailed analysis of the A2256 ROSAT PSPC data has been carried out by Henry, 
Briel, & Nulsen (1993). From their data, an isothermal model of the intra-cluster medium has 
been derived: 9 C = 5'.33 ± 0'.20, = 0.795 ± 0.020, n = 3.55 ± 0.18 x lO" 3 h 1 ' 2 cm" 3 . They derived 
a low-energy temperature from the ROSAT data of kT e = 6.9 ± 0.6 keV, which is consistent 
with the Ginga temperature. We choose to adopt the Henry et al. parameters and the David et 
al. Ginga temperature. 

6.1.4- Coma 

The Coma cluster is also a member of our sample. Briel, Henry & Bohringer (1992) 
present an isothermal X-ray model for the Coma gas with (5 = 0.75 ± 0.03 and 6 C = 10^5 ± 0f6 
(r c = 0.207 ± 0.012 h^ 1 Mpc for an assumed redshift of z = 0.0235). They adopt an electron 
temperature of kT e ff = 8.2 ± 0.2 keV, obtained from Ginga measurements. 

Earlier observations by Hughes, Gorenstein & Fabricant (1988) with EXOSAT give a slightly 
higher temperature of kT e ff = 8.5 ± 0.3 keV. After deprojection and subtraction of galactic 
absorption, they fit a model with (3 = 0.63 ± 0.03 and 9 C = 7'.Q ± 0f4, and central electron density 
of n w 3 x 10~ 3 /i 1/2 cm" 3 . For a best-fit model using the EXOSAT and TENMA data, they 
assume an isothermal core with the high temperature of kTi so = 9.1 ± 0.4 keV, which extends out 
to a radius of 6i so = 23' _ / , beyond which the temperature falls almost adiabatically (polytropic 
index 7 ~ 1.555). This gives a temperature profile outside the isothermal radius Ri so of 

-¥(7-D 

R > R lso . (39) 

This is the model adopted by Herbig et al. (1995), although with no, r c and (3 as given by Briel, 
Henry & Bohringer. 

The higher temperature may in fact be indicated by A SCA observations of Coma (see Fabian 
et al. 1994) which prefer temperatures of kT e ff ~ 9keV. For consistency with Herbig et al. we 
adopt the Briel, Henry & Bohringer no, r c and (3, and a temperature kT e ff = 9.1 ± 0.4keV. In 
the model calculations in this paper, we do not include an isothermal cutoff, and this makes only 
a few percent difference. 



T(R) = Tit 



1 + (R/r c ) 2 
1 + {R tso /r c ) 2 
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6.1.5. Model- Dependent Quantities 

Using the isothermal /3-model with electron temperature T e , the y at a point at projected 
radius t = R/r c from the cluster center has the familiar form 



y { t) = ^±. nQ a T r c f°° ds (l + t 2 + s 2 )'^ . = y (l + t 2 ) 1 '^ (40) 
ni e c J —oo ^ ' ^ ' 

where yo is the y-parameter at zero projected radius 

y = 7.12 X 10- ( J ) ( *U ( t*) . (41) 



r(f) 



10-3 h 1 / 2 cm" 3 J VIOkeVy yh^Mpc 



The dependence on h is due to the choice of the units for uq and r c , which in turn are determined 
from the X-ray measurements. The resulting factor of b}l 2 will be used to determine the value of 
the Hubble constant by comparison with the observed y-parameters. 

Given the model for the density profile n e (R,(j),Q in the cluster, we can determine the 
efficiency r] f JS at which our switched observations can recover the SZE that an ideal pencil-beam 
through the cluster center would measure 

Vobs = • (42) 

yo 

In addition, we can compute the efficiency rj g at which the SZE is measured with respect to an 
ideal Gaussian main beam 

USUI ,.„N 

Vg = — . (43) 
yg 

These efficiencies depend upon the model only through 9 C and (3. The derived efficiencies (for 
a pointing center at the cluster center R = 0) for the OVRO 5.5-m SZE observations are given 



in Table 1C. The uncertainties in the Ty f, s and r\ g were determined numerically using the stated 



uncertainties in the model 9 C and (3. 

Other model-derived quantities of interest are the equivalent spherical volume 

V s (R)=Aitrl [ R/r ° dtt 2 {l + t 2 )-^l 2 (44) 



and the Gaussian cylindrical volume 

rfi)rte) r oo r 2 t 2 

V g (L g ) = 2nrl V ^ J I dtteM-^)(l + t 2 ) 1/2 ^ /2 - (45) 

These quantities are the equivalent volumes for a uniform density cluster at the central density 
7100. These volumes are important for relating the observed X-ray emission and the observed SZE 
to the implied baryonic mass contained within the cluster. 



4 Note that equations ( fli] ) and ( fl5| ) can respectively be written in terms of the Incomplete Beta Function and 
the Incomplete Gamma Function. However, it is easiest to evaluate these integrals numerically, using Maple or 
Mathematica for example. 
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6.2. Baryonic Mass in Clusters 

Given knowledge about the electron temperature in the intracluster gas, we can use the y 
parameter to measure a baryonic mass M b for the ionized phase. For a general density model 



M b = J J J dRd(f)d(am b n e (R,(j),C). (46) 
Comparison with equation ( p4| ) gives for our cylindrical model 

TfL C 2 f 

M b = ^—a'm b 2ir / dRy(R), (47) 

where a' is the number of baryons per electron, and m b is the baryon (nucleon) mass. Considering 
only H and He at 12:1 in number of atoms, a' ~ 8/7. For a general temperature distribution, the 
effective temperature is given by 

rp _ U S dRd4>d{n e {RM)T e {R,<t>,Q 
eff ffJdRd<j>d(n e (R,<P,C) 1 ) 

which reduces to the single temperature T e ff = T e for an isothermal cluster medium at electron 
temperature T e . 



By combining (^7() with (pq), we can determine the temperature weighted baryonic mass 
within the cylinder defined by the 7^35 FWHM Gaussian beam of the 5.5-m telescope through the 
cluster center 



m e c 2 



M g = a'm b n V g = ^_ ^ 2irL g a'm b y g , (49) 



&TkT e ff 

where L g = D a 6 g as before. Hence, the baryonic mass may be written as 

M g = 4.407 x 10 14 ( ( - V ( h- 2 M . (50) 

9 \T eff J\h-iMpcJ \10- 5 Vg J ° v ; 

The efficiency factor rj g = y sw/llg converts the measured y sw into the y g within the Gaussian main 
beam. A better representation of the SZE in terms of a mass is the surface baryonic mass density 
within the Gaussian cylinder 



eff / V 1U 1 9 



The surface density is distance independent, and is a more consistent parameter than the mass, 
which will vary with the resolution L g . The values to use for y sw are the measured y-parameters 
ymeas-, with relativistic corrections applied using ([H]), found in Table 1C. The derived baryonic 
masses for our clusters are given in Table [ll]. The clusters have similar surface densities 
S s ~ 7 x lO 13 M Mpc- 2 . 

The factor r] g and the effective temperature T e ff are the only model-dependent quantities in 
M g and T, g . In Table 11, we use the X-ray model parameters listed in Table ||. Since ([5(]) and 
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(51) are linear in the observable y sw , the SZE is potentially a more accurate probe of the baryonic 



mass than the X-ray emission. 

The X-ray emission from clusters has been used to determine the baryonic mass fraction by 
comparison with derived total masses. This calculation has been done for Coma (White et al. 
1993) and A2256 (Henry, Briel, & Nulsen 1993). The SZE is an independent measure of the mass 
within the Gaussian cylinder of the beam. We can use the X-ray derived model to relate M g to 
the mass within the sphere 

M sze (R) = M g ^± (52) 

where V g is the Gaussian volume ( |45|) within the beam on the cluster. 

For our clusters A478, A2142, and A2256, as well as Coma, the baryonic masses are given 
in Table |ll]. In Table gravitational masses have been obtained from the literature and the 
baryonic fraction within some given fiducial radius Rq is computed. The model-dependent factors 
V s /Vg are listed for the assumed Ro, along with the 1-a uncertainties computed from the model 
uncertainties in r c and (5. 

White & Fabian (1995) discuss the "baryon overdensity" problem in the context of Einstein 
observations of a number of clusters, including A478 and A2142. Henry, Briel, & Nulsen (1993) 
give detailed models and masses for A2256, and White et al. (1993) compute the enclosed baryonic 
and gravitational masses for Coma. We discuss the results for each cluster, and the four clusters 
as a group, below. 



6.2.1. A478 

White & Fabian (1995) consider a radius of Rq = 0.976 h~ l Mpc within which they find an 
X-ray determined gas fraction M xray /Mt t = 0.091 ± 0.008 /i~ 3 / 2 . No uncertainties in the values for 
M to t are stated, although they are likely to be high (probably 20% or more). This should be kept 
in mind when evaluating the uncertainties for the clusters listed in the White & Fabian paper. 

For the model in Table ^, our SZE measurements give a Gaussian mass of 
M g = (2.58 ± 0.21) x 10 13 hr l M Q within the 5-m beam. This model gives a ratio 
V a /Vg = 2.99 ± 0.08 within R = 0.976 h' 1 Mpc, so we find an SZE-indicated baryonic 
mass of M sze = (7.71 ± 0.66) x 10 13 /i~ 2 Mq. Using the gravitational mass from White & Fabian, 
we get a baryonic fraction of M sze /Mt t = 0.166 db 0.014 h . 

The cluster A478 stands out with a higher baryonic fraction from both the X-ray and SZE 
measurements, and has a stronger SZE decrement than expected from the X-ray measurements 
(giving a lower implied Hubble constant from the ratio M xray /M sze ), when compared with the 
other clusters in this sample (see below, and in the next section). These discrepancies may be 
explained by elongation of the cluster along the line of sight, as indicated by its observed ellipticity 
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in the plane of the sky. We will discuss this further in the context of the Hubble constant in the 
next section. 



6.2.2. A2142 

For A2142, White & Fabian (1995) find an X-ray determined gas mass fraction 
M xray /M tot = 0.050 ± 0.003 h~^ 2 within R = 0.976 h' 1 Mpc. Our SZE measurements 
give a Gaussian mass of M g = (2.08 ± 0.22) x 10 13 h~ l M with an efficiency of 
V s /V y = 2.90 ± 0.43 for R = 0.976 h' 1 Mpc. Therefore, M sze = (6.03 ± 1.10) x 10 13 h~ 2 M Q and 
M sze /M tot = 0.060 ± 0.011 h- 1 . 



6.2.3. A 225 6 

Henry, Briel, & Nulsen (1993) fit a model with r c = 0.245 ± 0.009 ft," 1 Mpc, compared 
to the 5-m beam size of L g = 0.143 ft -1 Mpc at redshift z = 0.0581. They derive an enclosed 
mass of M tot = (5.1 ± 1.4) x 10 14 h~ l M within Rq = 0.76 h' 1 Mpc. From the X-ray data 
they find M xray /M tot = 0.063 ± 0.039 /i~ 3/2 . Assuming our isothermal model, we find a ratio 
Vs/Vg = 4.29 ± 0.05 within the sphere of radius Ro. The SZE measurements gave a Gaussian mass 
of M g = (7.1 ± 0.9) x 10 12 h- 1 M Q within the 5-m beam, thus M sze = (3.05 ± 0.39) x 10 13 h' 2 M Q 
in the sphere, and M sze /M tot = 0.060 ± 0.018 h' 1 . 



6.2.4. Coma 

For Coma, White et al. (1993) adopt a (model-dependent) total mass of 
Mtot = (1.10 ± 0.22) x 10 15 h' 1 M Q within a sphere of radius R = 1.5 h' 1 Mpc (the Abell radius). 
They find M xray /Mt t = 0.050 ± 0.013 /i~ 3 / 2 . Our adopted isothermal model with an assumed 
temperature of kT eff = 9.1 ±0.4keV gives the Gaussian mass of M g = (1.81 ±0.30) x 10 12 h' 1 M Q 
within the 5-m beam. For a spherical radius of Ro = 1.5 h^ 1 Mpc, the ratio V s /V g = 38.3 ±3.0, and 
thus M sze = (6.93 ± 1.27) x 10 13 hr 2 M & . Using this value, we find M sze /M tot = 0.063 ± 0.017 h' 1 , 
in agreement with the White et al. fraction for h = 0.62 ^[j' 3 g. 

For the Coma cluster, the 7f35 5-m beam is small (L g = 61.5 h^ 1 kpc) compared to the Abell 
radius Ro = 1.5 h -1 Mpc, and we are making a large Gaussian correction V s /V g ~ 38. In addition, 
the contribution of the SZE signal in the reference beams is significant, so the details of the 
electron temperature profile in these outer parts are more important than in the other clusters. 
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6.2.5. The Sample 

We see that in three of the four clusters (A2142, A2256, and Coma) the SZE determined 
baryonic fractions M sze /Mtot are consistent, with a mean of 0.061 ± 0.010 ft, -1 (the uncertainty 
from the individual error bars, not the scatter). We exclude A478 from this average due to the 
discrepancies between this cluster and the others in the sample (if A478 is included the mean 
becomes 0.087 ± 0.030 h~ l ). We should also include the 6.9% calibration uncertainty, giving 
(M sze /M tot ) = 0.061 ± 0.011 h-\ 

Strictly speaking, this is a lower limit on the baryon fraction Ms/Mtot, since we have 
not included the luminous mass in galaxies, and some of the dark matter may be baryonic. 
White et al. (1993) find a ratio M ga i/M tot = 0.009 ± 0.003 in Coma, compared to the fraction 
M xray /M tot = 0.050 ± 0.013 /i" 3/2 in hot gas. Henry, Briel, & Nulsen (1993) find similar relative 
fractions in galaxies and gas for A2256. Thus, we can safely assume that the luminous galaxies 
contribute around 20% or less of that mass contributed by the hot IGM. If we use the Coma value, 
and apply it to the sample as a whole, then {M B /M tot = 0.009 ± 0.003 + 0.061 ± 0.011 h~ l . Note 
that for the low values of h which are generally preferred, the contribution from luminous galaxies 
is further reduced relative to the gas. In this paper, we will use the SZE mass as a lower limit on 
the total baryon mass. 

Standard estimates of the fraction of closure density in baryons for homogeneous big-bang 
nucleosynthesis give 0.011 < S7#/i 2 < 0.015 (2d) (Smith, Kawano & Malaney 1993). However, 
recent measurements of the deuterium abundance in the Lyman-a forest clouds in several QSOs 
lead to incompatible values for Qb that lie outside this range. One group finds a low value for 
the deuterium abundance (Tytler, Fan and Buries 1996) implying a high fish 2 = 0.024 ± 0.006, 
while the other finds a high deuterium abundance (Rugers & Hogan 1996) implying a low 
tt B h 2 = 0.0062 ± 0.0008. 

We can now estimate the total mass density parameter 

_ Q.B < 

~ M B /M tot ~ M sze /M tot 

assuming that the baryonic fraction in clusters reflects that of the Universe as a whole. Using our 
reduced sample average (M sze /M tot ) = 0.061 ± 0.011 /r 1 , and assuming standard nucleosynthesis 
limits Q B h 2 = 0.013 ± 0.002, one obtains Q h < 0.21 ± 0.05 (and Q h < 0.15 ± 0.06 if A478 
is included in the mean). However, if one adopts the higher f2#/i 2 = 0.024 ± 0.006, then our 
data imply a significantly higher limit Qoh < 0.39 ± 0.12. On the other hand, if we use the 
high-deuterium value giving a low n B h 2 = 0.0062 ± 0.0008, then we find a low density parameter 

n h < o.io ±0.02. 

Using the standard nucleosynthesis values for we find the cluster data is consistent with 
Slo = 1 only for very low values of the Hubble constant (h ~ 0.21), or for similar values of baryonic 
mass segregation (^B,tot/^B,clus ~ 0.21). Neither of these is indicated by other cosmological data. 



(53) 



-32- 



Extremely low values of h are not consistent with the estimates derived by comparison with the 
X-ray emission except in the case of A478 (see below). On the other hand, this result is consistent 
with large-scale structure studies which yield values in the range 0.2 <; £l h ^ 0.3 (Efstathiou, 
Sutherland & Maddox 1990, Efstathiou, Bond & White 1992, Peacock & Dodds 1994). This 
widespread "baryon overdensity" problem has been seen in many clusters (e.g. White & Fabian 
1995). 

However, if the higher Qb from a low deuterium abundance is correct, then the values of the 
Hubble constant implied by comparison with the X-ray data (see below) would allow a critical 
density for the Universe. Conversely, adoption of the high deuterium, low baryon density Qb 
would only exacerbate the baryon overdensity problem, forcing us to accept a low density universe. 

Departures from isothermality or the coexistence of multiple phases in the intracluster 
medium will introduce errors in our determination, as will model errors in the extrapolation to the 
spherical masses. In particular, the determined values for the total binding masses are uncertain. 
Measurements of total mass surface density from weak gravitational lensing would be particularly 
well-suited to this method, as the angular size of a typical CCD frame is similar to that of the 
5.5-m beam width. 

Note that our SZE (and X-ray) measurements only count the baryons in the hot IGM. The 
luminous baryonic matter in galaxies, and any non-luminous baryonic matter (such as in brown 
dwarfs, Jupiters, or compact objects) would be in addition to this estimate. Thus, we place a 
lower limit on the total baryonic mass, and therfore the baryonic fraction. The baryon overdensity 
problem would only get worse if there were substantial contributions from these other baryon 
resevoirs. 

The difference between the SZE-based and X-ray based estimates of the baryonic fraction is 
due to the different dependences on the Hubble constant, h~ l versus /i -3 / 2 respectively. If we 
compare the X-ray and SZE numbers given above for A2142, A2256 and Coma (thereby excluding 
A478), the average ratio is (M xray /M sze ) = 0.889 ± 0.099 /i" 1 / 2 , and therefore in agreement for 
h = 0.79l°ir- We ex P lore this 

in more detail in the next subsection. 



6.3. The Hubble Constant 

The key to the determination of Hq lies in the observation that the X-ray brightness and 
SZE decrement scale differently with temperature and density. Because the core radius r c is 
determined from the observed angular size of the cluster (r c oc /i -1 ), the central densities of the 
cluster gas determined from X-ray data are proportional to /i 1 / 2 . The angular diameter-distance 
relation introduces a factor of h~ l into the y-parameter for the dependence on r c (see equation 
EJ). Thus, 

Dpredi the estimate of the switched measurements of the Compton y S w derived from the 
X-ray model in Table ||] using (^) and (^), are proportional to /i -1 / 2 , and the observations of the 
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actual SZE y mea s can therefore be used to find Hq: 

h= fVpre±\\ (54) 
\ Vmeas / 

The SZE model predictions y pre d, the measured ymeas-, and the inferred Hubble constant values 
are given in Table |l3[ The uncertainties in y pre d were computed from the model parameter 
uncertainties. These parameters were assumed to vary independently, though in fact they are 
correlated from the X-ray fitting procedure (particularly 6 C and (3). A more direct approach, 
comparing the X-ray data and SZE data directly, as in Birkinshaw et al. (1991) and Birkinshaw 
& Hughes (1994), would be preferable. As it is, using the available information, the uncertainties 
quoted here are likely to be slightly inflated, as the parameter correlations will reduce the overall 
uncertainty somewhat. The asymmetrical 1-cr error bars on h are computed from the symmetrical 
1-cr uncertainties on /i 1 / 2 . 

In the discussions below, it is clear that when a detailed examination of nearby clusters is 
made, significant departures from the spherically symmetric, smooth, isothermal cluster "ideal" 
are seen. Improved X-ray models from ASCA and ROSAT are critical to the use of the SZE to 
determine Hq. 



6.3.1. A478 

Using this X-ray model based upon the Allen et al. (1993) ROSAT observations and the 
Ginga temperature, we derive a Hubble parameter h 1/2 = 0.57 ±0.14 or h = 0.32+°;^. The largest 
contribution to the uncertainty comes from y pre d (±24%) rather than the 5.5-m measurement 
Vmeas (±7%). Thus, the largest improvement to be made is in the X-ray model. 

The low value of h = 0.32 is similar to the value we obtained when calculating the baryonic 
masses and mass-fractions in the previous section. A478 appears to have a much stronger SZE 
decrement than one would expect from the X-ray model, as well as a higher X-ray luminosity 
than one would expect from the size and velocity dispersion. One possible explanation for this is 
that A478 is significantly elongated along the line of sight, by around a factor of two compared 
with its dimensions in the plane of the sky. This would bring the implied value of the Hubble 
constant in line with the other clusters. It may also be that A478 is contaminated by the cooling 
flow emission. However, we have done some preliminary tests using a two-component model 
incorporating a low-temperature high-density phase, which gives nearly the same predicted SZE 
decrement, thus nearly the same Hubble constant. This cluster remains a puzzle and merits more 
detailed examination. 
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6.3.2. A2142 

The X-ray model predictions and observed SZE give = 0.69 ±0.26, or h = 0.48+°;^. The 
statistical error is dominated by the the uncertainty in y pre d (±43%) rather than in y meas (±6%). 
The largest single uncertainty is in the value of (3. It will be important to improve the model with 
a detailed analysis of the ROSAT data. 

6.3.3. A 225 6 

Using the Henry et al. parameters and the Ginga temperature, combined with the 5.5-m 
measurements of the SZE we derive h 1 / 2 = 0.85 ± 0.12, and thus h = 0.72+°; 22 >. The contributions 
to the error bar are ±7% from y pre d and ±12% from y meas - For this cluster, the dominant 
uncertainty is from the SZE measurement. It will be difficult to improve these measurements 
significantly, as it already has over 300 hours of integration time devoted to it. In addition, the 
CMB anisotropics on these scales are expected to be in the range 14 - 54 /iK rms (see § |j). Clusters 
with SZE decrements weaker than that in A2256 will be very difficult to use for determination of 
the Hubble constant. 

This cluster appears to have the best model, and cleanest X-ray and SZE data, though 
it is weaker than the others. Some possible problems not apparent in this analysis may be 
caused by the presence of the merging sub-clusters in the core, and the probable presence of 
very high-temperature gas indicated by preliminary reports from ASCA (J. P. Henry, private 
communication). This should be kept in mind when evaluating the A2256 data (and similarly for 
the other clusters), and in the not too distant future ASCA and ROSAT should be able to provide 
much better constraints on the cluster models. 

6.3.4. Coma 

We will adopt the Herbig et al. determination of the Hubble constant, rather than use our 
own isothermal model (see Table ^) . Because the Coma cluster subtends a large angle on the sky 
compared with the switching angle, there is significant contribution of the SZ decrement to the 
reference beams, and thus the details of where the isothermal cluster atmosphere cuts off makes a 
noticeable difference to the derived Hubble constant. 

Herbig et al. find a value of /i 1 / 2 = 0.843 ± 0.163, which after application of the 
relativistic correction (xrei = 1-029) and removal of the 6.9% calibration uncertainty, becomes 
h 1 / 2 = 0.819 ± 0.148. This gives us h = 0.67^Q22- Note that adoption of a lower temperature, 
such as the Briel et al. value, will reduce the derived Hubble constant for Coma. 
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6. 3. 5. Results for the Sample 

If the clusters in our sample are significantly ellipsoidal in shape, a value for Hq can only 
be obtained by averaging a number of results from the unbiased sample. We combine the 
measurements for A478, A2142, A2256, and Coma (using our relativistically corrected Herbig 
et al. value), and the average is listed in Table [13[ The most natural variable to average is 
h 1 / 2 = y pre d/y m eas: for which the measurement and model errors should enter, as nearly as 
possible, in a Gaussian fashion, and for which projection effects (see below) should average 
out in an orientation unbiased sample. In this case, we get the mean h 1 ! 2 = 0.733 ± 0.076, or 
Hq = 54^^ kms -1 Mpc -1 . Note that the reduced x 2 = 0.90 on the three degrees of freedom 
against the mean h 1 / 2 , so the spread in Hq is consistent with the (large) error bars. 

Up until this point, we have dealt with the "statistical" uncertainties introduced by the 
observations, calibration, and models. We should therefore now include the overall systematic 
calibration uncertainty of 6.9% (see §3.3). Because all of the observations were calibrated using 



the same scale (including Coma), any error is correlated between the four cluster measurements, 
and should thus be applied to the sample mean as a whole. Adding this uncertainty in quadrature, 
we find a sample average h 1 / 2 = 0.733 ± 0.091, or Hq = 54^^ kms -1 Mpc -1 . This is the value 
that we adopt. The reader is reminded also that if the individual cluster measurements are to be 
used from Table [ll] or from Table [l^, then the 6.9% calibration uncertainty should be added to 
the statistical error bars listed there. 

Because Hq depends upon the squares of the AT of the model and of the measurement, the 
fractional errors in each are effectively doubled before adding in quadrature to make up the error 
budget in the Hubble constant. Accurate determination of Hq therefore relies both upon accurate 
measurements of the SZE, and upon an accurate model of the state of the intracluster medium 
(see discussion in Birkinshaw et al. 1991, and Inagaki, Suginohara & Suto 1995). 

At the beginning of this section, we discussed the systematic errors introduced by the 
relativistic corrections to the SZE, and the effect of a cluster peculiar velocity. In the former case, 
corrections to the y-parameter were made, and in the latter case, the corrections were dismissed 
as unlikely to be important. 

The most serious potential source of systematic error in the determination of y pre d given the 
X-ray measurement is elongation of the cluster. Our analysis assumes that the line of sight extent 
of the cluster is the same as that in the plane of the sky. We have also assumed a spherical density 
profile in our analysis, though we would get the same result for an ellipsoidal cluster with r c as the 
geometric mean core radius. 

Deviations from circular symmetry in isophotes are not unusual: McMillan et al. (1988) 
studied 49 clusters observed by the Einstein satellite and found that the X-ray images had 
ellipticities of up to 0.5. As discussed in § [0|, the X-ray isophotes of A478 and A2256 show 
evidence for significant ellipticity. Elliptical isophotes on the sky imply, at least statistically, a 
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non-unity axial ratio in the line-of-sight dimension also. 

There is some indication that the cluster A478, with its high implied baryonic mass 
fraction and low implied Hubble constant, may be an example of a highly elongated cluster. 
If A478 were excluded from our sample average, then we would find h 1 / 2 = 0.788 ± 0.082, or 
Hq = 62^2 kms -1 Mpc" 1 . However, without any clear indication of a problem with A478 given 
the large error bars, we choose to adopt the entire sample average. 

It has been found that many clusters have significant cooling cores, and simple /3-models 
may be inadequate to describe the state of the gas in these cases. Edge and Stewart (1991a) list 
A478 and A2142 as cooling-flow clusters. Cooling-flow clusters are characterized by cores with 
gas at a significantly higher density and lower temperature than the surrounding gas. This causes 
a pronounced central peak in the X-ray surface brightness of the clusters. The cooling cores 
typically have radii of 50-200 kpc and temperatures reduced by up to a factor of four compared 
to the overall temperatures. Because the SZ decrement is more sensitive to the outer, low density 
regions of the gas distribution than is the X-ray flux, central densities calculated from Sx(r) at 
relatively large radii (Abramopoulos & Ku 1983, Jones & Forman 1984) should be used when 
possible. These densities are typically 2-3 times lower than the cooling core densities found by 
Edge and Stewart (1991a). Given X-ray images with high resolution and sensitivity, a better 
method would be to model n e and T e at large and small radii separately. Better models for the 
gas distribution are needed to account for the presence of cooling flows or other departures from a 
single spherically symmmetric smooth isothermal profile. 



7. Concluding Remarks 

In summary, we find significant detections of the SZE for the clusters A478, A2142 and A2256 
using the 5.5-m telescope at OVRO. These are the first observations to detect the effect in these 
clusters, although A478 and A2142 have been searched before (Lake & Partridge 1980, Birkinshaw 
et al. 1981, Birkinshaw & Gull 1984, Chase et al. 1987). Observations of contaminating radio 
sources were carried out on the OVRO 40-m telescope. When the SZE measurements of the X-ray 
flux-limited sample are complete, an orientation unbiased sample of clusters will be available for 
measurement of the Hubble constant. 

The SZE is a measure of the electron pressure in the ionized cluster medium — with 
knowledge of the electron temperature, the SZE is proportional to the baryonic mass in the IGM 
contained within the telescope beam. We find similar baryonic mass surface densities for the three 
clusters and Coma: £ 9 ~ 7 x 1O 13 M Mpc -2 . For A2142, A2256 and Coma, consistent estimates 
of the baryonic mass fraction M sze /M to t ~ 0.061 ± 0.010 hT 1 . This is a lower limit on the total 
baryonic mass, as the galaxies contribute M sze /Mt t ~ 0.009 ± 0.003. When compared with the 
standard primordial nucleosynthesis estimates for f2#, we find consistency between the SZE data 
and nucleosynthesis for cosmological density parameters in the range £loh ~ 0.21 ± 0.05. This 
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agrees with the values determined independently from large-scale structure and galaxy counts. 
The cluster A478 gives a discrepant (high) fraction of M sze /M to t ~ 0.17 /i -1 , and is likely elongated 
or heavily contaminated by the cooling flow (or both). 

Recent determinations of $1b using the deuterium abundances in Lyman-a absorption systems 
along the line-of-sight to QSOs give discrepant values higher and lower than the standard. If we 
adopt a high baryon density fish 2 = 0.024 ± 0.006 (Tytler, Fan and Buries 1996), then our data 
imply Q h < 0.39 ± 0.12. On the other hand, if Q B h 2 = 0.0062 ± 0.0008 (Rugers & Hogan 1996), 
then n h < 0.10 ± 0.02. 

By combining the measured SZE decrements with published X-ray models, we have determined 
the value of the Hubble constant implied for each of these clusters. Clusters A478 and A2142 
are complicated by excess core emission attributed to cooling flows, and better X-ray models 
for the gas distribution must be obtained. A2256 appears to be undergoing a merging event, 
and additional modeling must also be done here. With the preliminary models gleaned from the 
literature, we obtain an average of Hq = 54+^ kms" 1 Mpc" 1 for A478, A2142, A2256, and Coma. 
This average value tends toward the low side of the commonly accepted range, as do most of 
the other SZE determined values (e.g., Birkinshaw et al. 1991, Jones et al. 1993, Birkinshaw & 
Hughes 1994), though the large error bars place our measurement within 3d of practically all of 
the other values for Hq. 

Note that the high baryon density of Tytler, Fan and Buries (1996), our measurement of 
the baryon fraction, and our average value for the Hubble constant would imply a high overall 
density for the Universe 0,q < 0.72 ± 0.29. This is marginally consistent with a Universe with the 
critical density £1$ = \. However, adoption of the lower values of the baryon density would favor 
the acceptance of a low density parameter. 

Possible problems with using the SZE and X-ray measurements for inferring Hq include 
substructure in cluster atmospheres, elongation of clusters along the line of sight, and the presence 
of cooling cores. The SZE decrement is sensitive to the outer regions of the intracluster gas, which 
have not been well studied because of their relatively faint X-ray emission. Better models of the 
cluster atmospheres will soon be provided by the new generation of X-ray satellites, such as ASCA 
and AXAF. 

The advantage of determining Hq from a well-selected sample of clusters is the ability to use 
the distribution of derived Hq to test for variations in the astrophysical parameters of the cluster 
models assumed in the analysis. As this stage, our models are too uncertain and therefore our 
error bars too large, to assess any but the grossest deviations in derived h 1 / 2 . The most discrepant 
value is Hq = 32 ^ kms -1 Mpc -1 for A478, and even this is less than 2a from the mean. If 
massive clusters are inordinately elongated, estimates of Hq from individual clusters may be off 
by a factor of two or more. Exclusion of A478 from the sample average raises the value of Hq 
somewhat, though this step is unwarranted by the data at hand. 

The ability to recognize deviant clusters like A478 demonstrates the power of using this 
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sample of clusters. Clearly, completing the entire sample is the proper way to proceed using this 
method. However, the clusters reported here were selected as the first targets because they were 
free from strong contamination by radio sources. It will be very difficult to measure the SZE in 
the remaining clusters with the accuracy that we have been able to achieve with these first results. 

In the end, it will be the distribution of the Hq values for the sample that will tell us whether 
clusters are suitable tools with which to measure the expansion of Universe, or whether variations 
in shape and orientation, or density and temperature substructure introduce severe limitations in 
the determination of Hq by this method. Individual clusters are insufficient to make the case for 
one value of Hq or another, and it remains to be demonstrated that this method will yield reliable 
results. In any event, sounding of the intragalactic medium through combined SZE and X-ray 
measurements promises to provide important constraints upon multi-phase structures in the hot 
cluster atmosphere. This is as important as determining Hq, in that a number of cosmological 
tests, such as the baryon fraction Qb, rely upon observations of clusters of galaxies. 
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Table 1. The X-ray Selected Sample 



Cluster 


Position (B1950) 


z 


f x /w- u 


Wio 44 


kT 


@core 




R.A. 


Dec. 




erg cm~ 2 s -1 


h~ 2 ergs -1 


(keV) 


(arcmin) 


(Notes) 


[1] 


[1] 


[2] 


[2] 


[2] 


[3] 


[4] 


A85 


00:39:19.5 


-09:34:23 


0.0518 


6.37 


1.88 


6.2 


2f51 


A399 


02:55:07.6 


+12:50:47 


0.0715 


3.41 


1.94 


5.8 


I'm 


A401 


02:56:12.0 


+13:22:43 


0.0748 


5.88 


3.68 


7.8 


4^68 


A478 


04:10:40.1 


+ 10:20:21 


0.0900 


6.63 


6.02 


6.6 


lf85 


A754 


09:06:49.7 


-09:28:57 


0.0528 


8.53 


2.62 


9.1 


8^45 


A1651 


12:56:48 


-03:55:00 


0.0825 


3.67 


2.80 


7.0 




Coma 


12:57:19 


+28:13:24 


0.0232 


25.4 


1.49 


8.1 


10f5 


A1795 


13:46:35.4 


+26:50:23 


0.0616 


5.30 


2.23 


5.3 


3^03 


A2029 


15:08:27.2 


+05:55:56 


0.0767 


7.52 


4.92 


7.8 


lf58 


A2142 


15:56:15.8 


+27:22:38 


0.0899 


7.50 


6.80 


8.7 


3^69 


A2256 


17:06:56.3 


+78:43:02 


0.0601 


5.20 


2.08 


7.5 


5^33 



Note. — [1] Positions from Einstein IPC, except A1651 from Abell, et al. (1989). [2] Redshifts, 
X-ray fluxes and luminosities (2-10 keV) from Edge et al. (1990). Luminosities assume qo = 1/2. 
[3] X-ray temperatures from Einstein MPC (David et al. 1993) except: A478 from Ginga and 
ROSAT (Allen et al. 1993); A1651 from HEAO-1 (Edge et al. 1990); Coma, A1795, A2142, 
A2256 from Ginga (David et al. 1993). [4] Core radii from Jones & Forman (1984) except: A401, 
A754, A2142 from Abramopoulos & Ku (1983); A478 from Allen et al. (1993); Coma from Briel 
et al. (1992); A2256 from Henry et al. (1993). 



Table 2. 



Pointing Positions for SZE Observations 



Cluster Field Position (J2000) 





R.A. 


Dec. 


A478L 


03:54:25.01 


+10:27:40.73 


A478 


04:13:25.01 


+10:27:40.73 


A478T 


04:32:25.01 


+10:27:40.73 


A2142L 


15:43:18.00 


+27:13:32.00 


A2142 


15:58:18.00 


+27:13:32.00 


A2142T 


16:13:18.00 


+27:13:32.00 


A2256L 


16:48:54.50 


+78:38:27.00 


A2256 


17:03:54.50 


+78:38:27.00 


A2256T 


17:18:54.50 


+78:38:27.00 



Table 3. OVRO 5.5-m Standard Calibration Scale 



Scale 


TMars 


SDR21 


•SV7027 


53C286 


Relative to Jupiter 


1.252+0.022 


134.7+2.1 mJy/K 


35.5+1.1 mJy/K 


13.2+0.6 mJy/K 


T j upi ter = 144 ± 8 K 


180+10 K 


19.4+1.1 Jy 


5.11+0.33 Jy 


1.90+0.14 Jy 
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Table 4. 


OVRO 5.5-m 


Measurements of SZE 




A478 


A2142 


A2256 




( M K) 


M K) 


( M K) 


M-(L+T)/2 


-375 ± 24 


-420 ± 19 


-218 ± 14 


M-L 


-366 ± 27 


-399 ± 20 


-217 ± 15 


M-T 


-386 ± 27 


-451 ± 22 


-238 ± 16 


MAIN 


-139 ±20 


-214 ± 17 


-310 ± 12 


LEAD 


244 ± 21 


159 ± 17 


-94 ± 12 


TRAIL 


242 ± 20 


226 ± 16 


-82 ± 12 


L-T 


-13 ±27 


-66 ± 20 


-22 ± 15 


(L+T)/2 


248 ± 15 


188 ± 13 


-74 ± 13 


v eff 


800 


840 


2020 




1146 


1337 


4117 


t (hrs) 


64 


85 


310 
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Table 5. Sources within 9' of field centers or reference arcs. 



Source 


Position (J2000) 




Slow 


-518.5 


a 


5*32 




K.A. 


Dec. 




(mJy) a 


(mJy) 




(mJy) 


A478L.1 


03:54:57.2 


+ 10:12: 


:30 


33+7 


< 3.9 


< -1.4 


< 1.8 


A 478 1 


04:13:34 


+ 10:28: 


■04 


1 ft_|_7 b 


<r Q 1 






A2142L.1 


15:41:46.8 


+27:05: 


:54 


51+8 


13.5+1.1 


-0.99 + 0.13 


7.8+0.9 


A2142L.2 


15:42:58.4 


+27:06: 


:46 


34+7 


< 5.1 


< -1.2 


< 2.6 


A2142.1 


15:57:11.2 


+26:51: 


:31 


52+9 


8.5+1.9 


-1.35 + 0.21 


4.1+1.1 


A2142.2 


15:58:14.3 


+27:15: 


:48 


44+8 


7.9+2.0 


-1.28 + 0.23 


3.9+1.1 


A2142.3 


15:58:47 


+27:18: 


:06 


18+7 c 


< 6.3 


< -0.3 


< 5.4 


A2142.4 


15:59:05 


+27:03: 


:19 


42+15 c 


< 5.1 


< -0.9 


< 3.2 


A2142T.1 


16:12:26.3 


+27:23: 


:16 


49+8 


14.3+1.5 


-0.92 + 0.14 


8.6+1.1 


A2256.1 


17:02:09 


+78:40: 


:56 


48+8 d 


< 6.9 


< -0.7 


< 4.8 


A2256.2 


17:02:28 


+78:42: 


:57 


166+11 d 


< 4.8 


< -1.3 


< 2.3 


A2256.3 


17:03:03 


+78:36: 


:40 


62+4 d 


< 5.1 


< -0.9 


< 3.0 


A2256.4 


17:03:09 


+78:40: 


:00 


39+3 d 


< 6.0 


< -0.7 


< 4.1 


A2256.5 


17:03:28 


+78:37: 


:58 


157+10 d 


8.3+2.5 


-1.14 + 0.12 


4.4+1.4 


A2256.6 


17:03:51 


+78:46: 


:03 


185+13 d 


< 10.8 


< -1.1 


< 6.0 


A2256.7 


17:04:48.9 


+78:38: 


:29 


11+1 d 


< 10.5 


< 0.03 


< 10.7 



a Flux density at 4.85 GHz from 87GB unless otherwise noted. 
b Flux density at 10.7 GHz (Andernach et a/.1986). 
c Flux density at 2.7 GHz (Andernach et a/.1986). 
d Flux density at 1.415 GHz (Bridle et oi.1979). 
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Table 6. Corrections for sources within 9' of field centers or reference arcs. 



Source 


Rad 


Pos 
Ang 


d 




( uK) 


AT 

z -* J - max 

( uK) 


A478L.1 


17:i2 


152° 


5:81 


62° 


< 11 


< 1 


A478.1 


2:24 


80° 






< 57 


< 44 


A2142L.1 


21f67 


249° 


0:49 


-21° 


48±6 


-24 ±3 


A2142L.2 


8f05 


213° 






< 16 


< 1 


A2142.1 


26?57 


214° 


4:41 


-56° 


26±7 


-5±1 


A2142.2 


2f41 


340° 






24±7 


18 ±5 


A2142.3 


7:90 


55° 






< 34 


< 1 


A2142.4 


14^62 


134° 


7:54 


44° 


< 20 


< 1 


A2142T.1 


15^05 


310° 


7:il 


40° 


53±7 


-2.0 ±0.3 


A2256.1 


5:75 


296° 






< 30 


< 5 


A2256.2 


6:i9 


317° 






< 14 


< 2 


A2256.3 


3:i0 


235° 






< 19 


< 11 


A2256.4 


2:72 


305° 






< 26 


< 17 


A2256.5 


i:39 


250° 






27±9 


25 ±8 


A2256.6 


7:60 


359° 






< 37 


< 2 


A2256.7 


2:68 


89° 






< 67 


< 46 



Note. - - Radius (arcmin) and Position Angle are measured 
from center of main field. Distance d (arcmin) is measured from 
center of reference arc at parallactic angle ip p . 
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Table 7. Total source contributions to 5.5-m SZE Fields 



Field 


N on 


AT on 


Nref 


max AT re f 








(MK) 




(fiK) 


(/xK) 


A2142L 


1 


0.6 


1 


-23.7 


-0.7 


A2142 


1 


17.8 


1 


-4.8 


16.2 


A2142T 





0.0 


1 


-2.0 


-0.2 


A2256 


1 


24.5 





0.0 


24.5 



Table 8. Final SZE Results 



Cluster 


AT 5 . 5m (obs) 


AT 5 . 5m (src) 


AT 5 . 5m (corr) 




(fjK) 


( M K) 


(/xK) 


A478 


-375 ± 24 


<44 


-375 ± 28 a 


A2142 


-420 ± 19 


17 


-437 ± 25 


A2256 


-218 ±14 


25 


-243 ± 29 



a A478 source limit 3a, used la for uncertainty. 



-47- 



Table 9. X-Ray Cluster Parameters 



Cluster 


z 


Score 


kT e 


no/10- 3 


P 








(arcmin) 


(keV) 


(h 1 ' 2 cm- 3 ) 




(h- 1 Mpc) a 


A478 


0.0881 


1.93±0.30 


6.56±0.09 


9.55±1.75 


0.667±0.029 


0.128±0.020 


A2142 


0.0899 


3.69±0.14 


8.68±0.12 


6.97±0.41 


1.0±0.3 


0.249±0.009 


A2256 


0.0581 


5.33±0.20 


7.51±0.11 


3.55±0.18 


0.795±0.020 


0.245±0.009 


Coma 


0.0235 


10.50±0.60 


9.10±0.40 


4.09±0.06 


0.750±0.030 


0.207±0.012 



a Assumes qo = 1/2. 



Note. - X-ray temperatures from Ginga (David et al. 1993), except Coma (Hughes, 
Gorenstein & Fabricant 1988). Other parameters: A478 from Allen et al. (1993) and our 
own analysis of ROSAT data; A2142 from Abramopoulos & Ku (1983); A2256 from Henry 
et al. (1993); Coma r c , uq and (5 from Briel, Henry & Bohringer (1992). 



Table 10. Efficiencies and Relativistic Corrections 



Cluster 


Vg 


Vobs 


Xrel 


Umeas 










(io- 5 ) 


A478 


0.826±0.028 


0.414±0.008 


1.022 


7.41±0.55 


A2142 


0.947±0.081 


0.498±0.090 


1.030 


8.70±0.50 


A2256 


0.818±0.014 


0.612±0.003 


1.026 


4.82±0.58 


Coma 


0.618±0.026 


0.563±0.021 


1.029 


6.13±0.93 



-48- 



Table 11. Baryonic Masses from the SZE 



Cluster 




S 9 


Mg 




(h- 1 Mpc) a 


(10 13 M o Mp C - 2 ) 


(10 13 h- 2 M ) 


A478 


0.207 


9.59 ±0.79 


2.58 ±0.21 


A2142 


0.211 


7.42 ± 0.77 


2.08 ±0.22 


A2256 


0.143 


5.50 ±0.67 


0.71 ±0.09 


Coma 


0.0615 


7.64 ± 1.25 


0.18 ±0.03 



a Assumes qo = 1/2. 



Table 12. Baryonic Fractions in the Clusters 



Cluster 


i?o 


v,m/v g 


M sze 


M tot 


M sze /M tot 




(h^ 1 Mpc) 




(W 13 h- 2 M Q ) 


(10 13 h- l M & ) 




A478 a 


0.976 


2.99±0.08 


7.7 ±0.7 


46.4 


0.166 ±0.014 


A2142 a 


0.966 


2.90±0.43 


6.0 ± 1.1 


100.5 


0.060 ±0.011 


A2256 b 


0.76 


4.29±0.05 


3.0 ±0.4 


51 ± 14 


0.060 ± 0.018 


Coma c 


1.50 


38.3±3.0 


6.9 ± 1.3 


110 ±22 


0.063 ± 0.017 



a A478 and A2142 M tot from White & Fabian (1995). No uncertainties given. 
b A2256 M tot from Henry Briel k Nulsen (1993). 
c Coma M tot from White el al. (1993). 
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Table 13. Hubble Constant from the SZE 



Cluster 


Vmeas 

(10- 5 ) 


Upred 

(io- 5 h- 1 ' 2 ) 




H 

(km s _1 Mpc -1 ) 


A478 
A2142 
A2256 
Coma 


7.4 ±0.6 

8.7 ±0.5 

4.8 ±0.6 
6.1 ±0.9 


4.2 ± 1.0 

6.0 ±2.2 

4.1 ±0.3 


0.57 ±0.14 
0.69 ± 0.26 
0.85 ±0.12 
0.82 ±0.15 


on+18 

48 + 43 

^°-29 
79 +22 
' Z -19 

u ' -22 


Sample 






0.73 ± 0.08 


54 +12 
d -ii 
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Fig. 1. — The effect of the data editing and filtering method upon the final referenced results. We 
show results for twenty different editing methods, for both unweighted (dashed lines) and weighted 
(solid lines) averages. The points are plotted slightly displaced from the filter method. The methods 
are roughly in order of fraction of data accepted, ranging from 100% (method 1) to 44% (method 20). 
The effective fraction of data used in the weighted averages ranges from 37% to 27% across the filter 
methods. We adopt method 15 with weighting (marked with the triangle), in which 57% of the data is 
retained, with an effective weighted fraction of 31%. 



Fig. 2. — Dependence of the LEAD, MAIN, and TRAIL (upper) and MLT and L-T referenced data 
(lower) on parallactic angle are shown for A478. In the upper panels, the LEAD fluxes are offset to the 
left of the proper ip p and the TRAIL fluxes are offset to the right. In the lower panel, the solid and 
dotted horizontal lines depict the means of the MLT and L-T. No source contributions were subtracted 
from A478. 



Fig. 3.— Dependence of the LEAD, MAIN, and TRAIL (upper) and MLT and L-T referenced data 
(lower) on parallactic angle are shown for A2142. In the lower panel, the solid and dotted horizontal 
curves represent response of the MLT and L-T data to the contaminating sources. Note the signal at 
■tpp = —21° that matches the deviation in the data. The expected signal from the source is below that 
observed, and we may have underestimated its contribution. However, because the number of data 
points in this ip p range is small, this does not affect the results significantly. 



Fig. 4. — Dependence of the LEAD, MAIN, and TRAIL (upper) and MLT and L-T referenced data 
(lower) on parallactic angle are shown for A2256. In the lower panel, the dotted horizontal line depicts 
the expected zero mean L-T, as no scans of these fields were available. The dotted line is the MLT 
mean before main beam source correction, and the solid line is the mean after source correction. 
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Fig. 1. — The dependence of SZE results on data editing. 
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Fig. 2.— SZE in A478. 



-53- 



Fig. 3.— SZE in A2142. 
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Fig. 4.— SZE in A2256. 



